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ABSTRACT 

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) was a multi-cycle treasury 
program on the Hubble Space Telescope ( HST) that surveyed a total area of ~0.25 deg 2 with ~900 HST orbits 
spread across five fields over three years. Within these survey images we discovered 65 supernovae (SNe) of all 
types, out to z ~ 2.5. We classify ~24 of these as Type la SNe (SNe la) based on host galaxy redshifts and SN 
photometry (supplemented by grism spectroscopy of six SNe). Here we present a measurement of the volumetric 
SN la rate as a function of redshift, reaching for the first time beyond z = 2 and putting new constraints on SN la 
progenitor models. Our highest redshift bin includes detections of SNe that exploded when the universe was only 
~3 Gyr old and near the peak of the cosmic star formation history. This gives the CANDELS high redshift sample 
unique leverage for evaluating the fraction of SNe la that explode promptly after formation (<500 Myr). Combining 
the CANDELS rates with all available SN la rate measurements in the literature we find that this prompt SN la 
fraction is /p=0.53 sta ^ 9 s ± 0 0 i 6 °, consistent with a delay time distribution that follows a simple t 1 power law for 
all times t > 40 Myr. However, mild tension is apparent between ground-based low-z surveys and space-based 
high-z surveys. In both CANDELS and the sister HST program CLASH (Cluster Lensing And Supernova Survey 
with Hubble ), we find a low rate of SNe la at z > 1. This could be a hint that prompt progenitors are in fact 
relatively rare, accounting for only 20% of all SN la explosions — though further analysis and larger samples will 
be needed to examine that suggestion. 

Key words: infrared: general - supernovae: general - surveys 
Online-only material: color figures 


1. INTRODUCTION 

The prevailing model for a Type la supernova (SN la) pro- 
genitor system begins with a binary system in which the pri- 
mary star evolves to become a white dwarf (WD). The WD 
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acquires mass from its companion star, approaches the Chan- 
drasekhar limit, and explodes in a thermonuclear runaway (see 
Hillebrandt & Niemeyer 2000; Livio 2001, for reviews). The 
companion star that feeds the WD and thereby sets off the ther- 
monuclear bomb is one of the key components of this model, but 
remains a topic of ongoing debate. In single degenerate (SD) 
models, the companion is a main sequence or evolved giant 
star, transferring mass via Roche lobe overflow, stellar winds, 
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Figure 1 . Volumetric SN la rates before completion of the CANDELS and 
CLASH SN surveys. Assorted ground-based surveys are plotted as white circles 
(Blanc et al. 2004; Botticella et al. 2008; Cappellaro et al. 1999; Dilday et al. 
2010; Hardin et al. 2000; Horesh et al. 2008; Graur & Maoz 2013; Li et al. 
2011; Melinder et al. 2012; Pain et al. 2002; Perrett et al. 2012; Rodney & 
Tonry 2010; Tonry et al. 2003). Three high redshift SN surveys are highlighted: 
gray circles for the Subaru Deep Field (SDF; Graur et al. 201 1), blue downward 
triangles for volumetric (not cluster) rates from the Cluster Supernova Survey 
(CSS; Barbary et al. 2012), and green upward triangles for the GOODS and 
PANS surveys (Dahlen et al. 2008). 

(A color version of this figure is available in the online journal.) 


or other means (Whelan & Iben 1973). In double degenerate 
(DD) models the companion is another WD, merging with the 
primary after a period of orbital decay driven by gravitational 
wave radiation (Iben & Tutukov 1984; Webbink 1984). More re- 
cent variations on these pathways to explosion include the “core 
degenerate” scenario (Kashi & Soker 2011) and perturbation- 
induced mergers in triple systems (Thompson 2011). 

The SN la explosion rate as a function of redshift, SNR(z), 
can provide an important observational test to constrain SN la 
progenitor models and possibly distinguish between them. In 
this paper we will present measurements of the SN la rate as a 
function of redshift and use them to place new constraints on 
SN la progenitor models, particularly on the fraction of SN la 
progenitors that explode within 500 Myr after their formation. 

Suppose we have a burst of star formation in a galaxy, such 
that the star formation rate can be approximated by a delta 
function in time. Binary population synthesis modeling gives us 
the initial conditions of all the binaries (mass, orbital separation, 
etc.), and a progenitor model sets the conditions necessary for 
explosion as an SN la. Using a stellar evolution model, one can 
follow the binary systems as they evolve, measuring the delay 
time distribution (DTD) between formation and explosion. To 
put constraints on SN la progenitor models, we can translate this 
DTD to cosmic scales and compare it to the observed volumetric 
SN la rate as a function of look-back time, as first proposed by 
Madau et al. (1998). 

As shown in Figure 1 , recent measurements of the SN la rate at 
low redshift (z < 1) are in good agreement, consistently finding 
that the SNR(z) rises steadily to at least z ~ 1 (e.g., Rodney & 


Tonry 2010; Dilday et al. 2010; Perrett et al. 2012). However, 
at z > 1 the trend of the SNR(z) curve is much less clear. The 
spectral energy distribution of an SN la peaks in the rest-frame 
B band with an absolute magnitude around —19.5. At z = 1.2 
that peak brightness becomes fainter than 25th magnitude in the 
observer’s z band — making discovery and light curve follow-up 
nearly impossible for ground-based observatories. 

For that reason, space-based surveys using the Hubble Space 
Telescope's Advanced Camera for Surveys (HST+ACS) have 
been the primary vehicle for tracking the SNR(z) to z ~ 1.5. 
The GOODS+PANS surveys were the first programs to extend 
rate measurements beyond z ~ 1 (Dahlen et al. 2004; Dahlen 
et al. 2008), and their measured rates suggested a peak in the SN 
la rate at z ~ 1.2, with a decline at higher redshifts. Independent 
examination of the same survey data recovered the same trend 
(Kuznetsova et al. 2008), although both analyses were limited by 
a small sample size in the highest redshift bin. Subsequently, the 
Cluster Supernova Survey (CSS) of the Supernova Cosmology 
Project used HST+ACS to measure the volumetric SN la rate 
(Barbary et al. 2012). These data revealed a similar peak and 
decline, although with even larger uncertainty in the high-z 
bins. From the ground, the Subaru Deep Field (SDF) SN survey 
used the Suprime-cam imager on the Subaru telescope to reach 
similar redshifts (Poznanski et al. 2007; Graur et al. 2011). As 
can be seen in Figure 1, these SDF rates formally show no 
decline in the highest redshift bin, but they are consistent with 
the HST+ACS results, within the errors. 

The HST-ACS high-z SNe la generally have reliable clas- 
sifications, based on well-sampled multi-band light curves, 
spectroscopic redshifts, and HST grism spectroscopy of most 
SN la candidates. However, due to the relatively small survey 
area, these programs have very large statistical uncertainties 
(Dahlen et al. (2008) have SNe in their highest redshift bin; 
Barbary et al. (2012) have ~1). In contrast, the SDF survey 
built up a larger sample (10 SNe la at z ~ 1.5) but their sur- 
vey design introduced potential for large systematic biases. The 
SDF epochs were spaced by ~1 year, meaning that the phase 
of the SN light curve at discovery was unconstrained, and the 
classification of detected SNe was based on only a single epoch 
of photometric data in the R, i', z! bands. Furthermore, red- 
shifts for the SDF high-z SN sample were based almost exclu- 
sively on photometric redshift estimates of the SN host galax- 
ies, not as precise or reliable as spectroscopic redshifts — but see 
Frederiksen et al. (2014) for one spectroscopic confirmation of 
an SDF host galaxy at z = 1 .55. 

An apparent peak in the SN la rate at z ~ 1 and a 
decline toward z = 1.5 has been interpreted as indicating a 
delay of > 1 Gyr between formation and explosion for most 
SN la (Strolger et al. 2004; Strolger et al. 2010). This would 
be broadly consistent with some SD models, and inconsistent 
with DD models, which typically predict a large fraction of SNe 
la that explode promptly after star formation (within 1 Gyr). A 
clear measurement of the shape of the SN la rate function at 
z > 1 would provide an important constraint on DTD models, 
and would go a long way toward resolving the question of 
whether an SD or DD model could be the dominant progenitor 
channel for all SNe la at all redshifts. Given the problems with 
current high-z SN rates, there is a clear need to improve the 
measurement by expanding the sample of well-classified SNe 
at z > 1. 

In this paper we present a measurement of the SNR(z) from 
a sample of 65 SNe discovered in the CANDELS SN program, 
extending the SNR(z) measurement for the first time to z = 2.5. 
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CANDELS SN Survey Fields a 
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Field 

R.A. 

(J2000) 

Decl. 

(J2000) 

WFC3-IR 

Tiles/Epoch 

Searchable Area 
(arcmin 2 ) 

SN Search 
Epochs (MJD) b 

COSMOS 

10:00:28 

+02:12:04 

44 

196.8 (9' x 22') 

[55905], 55953 

EGS-A 

14:19:18 

+52:49:30 

25 

106.2 (i of 7' x 32') c 

[55653], 55703 

EGS-B 

14:19:18 

+52:49:30 

20 

92.9 ( j of 7' x 32') c 

[56387], 56437 

UDS 

02:17:38 

-05:12:00 

44 

207.1 (9' x 22') 

[55512], 55562 

GOODS-S Wide 

03:32:42 

-27:53:37 

~8 

39.4 (4' x 10') 

[55573], 55621 

GOODS-S Deep 

03:32:28 

-27:46:01 

~15 

66.5 (7' x 10') d 

[55480], 55528, 55578, 55624, 55722, 
55774, 55821, 55860, 55921, 55974 

GOODS-N Wide NE 

12:37:29 

+62:18:40 

~8 

38.1 (4' x 10') 

[56183], 56238 

GOODS-N Wide SW 

12:36:20 

+62:10:25 

~10 

49.5 (5' x 10') 

[56020], 56073 

GOODS-N Deep 

12:36:55 

+62:14:19 

~15 

66.8 (7' x 10') d 

[56020], 56073, 56126, 56183, 56238, 
56297, 56348, 56402, 56458, 56511 


Notes. 

a Coordinates give approximate center of each CANDELS IR survey field. 
b Mean date of observation epoch. First epoch listed [in brackets] provided IR template images. 

c The CANDELS EGS field was divided into two interlocking halves, observed separately in 201 1 and 2013. See Grogin et al. (2011) for details. 
d The deep field search areas vary by epoch. The given value reflects the average. 


This SN survey is a joint operation of two HST Multi-Cycle 
Treasury (MCT) programs: the Cosmic Assembly Near-infrared 
Deep Extragalactic Legacy Survey (CANDELS; Pis: Faber and 
Ferguson; Grogin et al. 2011; Koekemoer et al. 2011), and the 
Cluster Lensing and Supernovae search with Hubble (CLASH; 
PI: Postman; Postman et al. 2012). The SN discovery and follow- 
up for both programs were allocated to the HST MCT SN 
program (PI: Riess). The results presented here are based on 
the full five fields and ~0.25 deg 2 of the CANDELS program, 
observed from 2010 to 2013. A companion paper presents the 
SN la rates from the CLASH sample (Graur et al. 2014). A 
composite analysis that combines the CANDELS+CLASH SN 
sample and revisits past HST surveys will be presented in a 
future paper. 

In Section 2 we describe the SN search component of the 
CANDELS survey, and in Section 3 we describe our detection 
efficiency measurements. Our photometric SN classifications 
are presented in Section 4, properties of the SN host galaxies 
are described in Section 5, and in Section 6 we detail new 
grism spectroscopy for four of our SNe. The rate calculation 
is described in Section 7 and we discuss the consequences for 
SN la progenitor models in Section 8. Finally, a summary is 
presented in Section 9. In tables and figures throughout the 
paper, we present the subset of 14 SNe with z > 1 .5 in the main 
body of the text, with the remaining 51 shown in Appendix B. 
Throughout this work we assume a flat ACDM cosmology with 
Ho ~ 70, Q m = 0.3, and Qa = 0.7. 

2. THE CANDELS SN SURVEY 

The three-year CANDELS program was designed to probe 
galaxy evolution out to z ^ 8 with deep infrared (IR) and optical 
imaging of five well-studied extragalactic fields: GOODS-S, 
GOODS-N, COSMOS, UDS, and EGS. 26 As described fully 
in Grogin et al. (2011), the CANDELS program includes both 
“wide” and “deep” fields. The wide component of CANDELS 
comprises the COSMOS, UDS, and EGS fields, plus one-third of 
the GOODS-S field and one half of the GOODS-N field — a total 

26 GOODS-S/N: the Great Observatories Origins Deep Survey South and 
North (Giavalisco et al. 2004); COSMOS: the Cosmic Evolution Survey 
(Scoville et al. 2007; Koekemoer et al. 2007); UDS: the UKIDSS Ultra Deep 
Survey (Lawrence et al. 2007; Cirasuolo et al. 2007); EGS: the Extended 
Groth Strip (Davis et al. 2007). 


Table 2 

Typical Exposures for a Single SN Search Epoch 


Camera 

Filter 

Exposures 
(Aexp x s) 

Limiting 

Magnitude a 

WFC3-IR 

F160W (H) 

2 x 600 

25.4 

WFC3-IR 

F125W (7) 

2 x 500 

25.8 

WFC3-UVIS 

F350LP (W) 

1 x 430 

27.8 

ACS-WFC 

F814W (7) 

2 x 700 

27.3 

ACS-WFC 

F606W (V) 

2 x 350 

28.1 


Note. a Vega magnitude that yields S/N ~ 5 in the given exposure sequence. 


survey area of 730 arcmin 2 . The CANDELS survey provides 
two visits to each wide field, spaced by ~50 days. The “deep” 
component of CANDELS came from the central 67 arcmin 2 
of each of the GOODS-S and GOODS-N fields. These deep 
regions were each visited 15 times over the course of two 
years (2010-2012 for GOODS-S, 2012-2013 for GOODS-N). 
Only 10 of those visits are used for SN discovery (the other 
visits lack template data for generating difference images), and 
those 10 epochs are also spaced at a cadence of ~50 days. The 
CANDELS fields analyzed in this work are described in Table 1 . 

Table 2 presents the exposure times and 5a limiting magni- 
tudes for a typical single-epoch set of exposures. Each CAN- 
DELS visit includes a set of four infrared exposures from the 
Wide Field Camera 3 (WFC3) IR detector: two in F160W 
(H band) and two in F125W (7 band). These are the search 
filters for the CANDELS SN survey (i.e., all SNe in our sample 
are infrared detections). Additionally, each observation set in- 
cludes a broad optical band, which helps to distinguish SNe la 
from Core Collapse Supernovae (CCSNe) and other transients 
(see Section 4). In ~80% of the SN search visits, this blue com- 
ponent is collected within minutes of the IR exposures as a single 
exposure using the WFC3 UVIS camera in the F350LP filter (a 
broad “white light” filter that we refer to as the “Wband”). In 
the remaining ~20% of visits (in the wide fields) the W band 
exposure is replaced with ACS observations in the F606W filter 
(broad V band), and complemented by the ACS F814W filter 
(broad /band). These ACS observations come from coordinated 
parallel visits and are taken within three days of the primary 
IR visit. 
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In addition to the ~750 HST orbits devoted to survey imag- 
ing in the CANDELS program, an additional 150 HST orbits 
were allocated for target of opportunity (ToO) follow-up obser- 
vations of newly discovered SNe. Another 52 orbits were pro- 
vided by the CLASH program, so the total CANDELS+CLASH 
SN follow-up allocation was 202 orbits. These follow-up visits 
provided supplementary imaging and slitless spectroscopy ob- 
servations to aid in the classification of SN candidates, and to 
measure the light curves of SNe la, allowing distance determi- 
nations for cosmology. 

2.7. Data Processing Pipeline 

All CANDELS survey images were processed through a data 
processing pipeline optimized for the detection of SNe by human 
searchers. This pipeline is similar in function to the CANDELS 
and CLASH pipelines (Koekemoer et al. 2011; Postman et al. 
2012), but includes some important differences specific to the 
SN search. There are four principal components in the pipeline: 
calibration, image combination, template subtraction, and fake 
SN planting. 

In the calibration stage, RAW images from HST are processed 
into FLT images using the STSDAS calibration tools provided 
by the Space Telescope Science Institute. 2 This includes bias 
correction, dark subtraction, flat fielding, and “up-the-ramp” 
fitting for cosmic ray rejection, as appropriate for each camera 
and detector. 

The image combination step uses the MultiDrizzle soft- 
ware (Koekemoer et al. 2002; Fruchter & Hook 2002) to com- 
bine multiple dithered images in the same filter from the same 
observing epoch, while also removing the geometric distortion 
of the HST focal plane. For each drizzled WFC3-IR image, we 
then generate a template image that combines all intersecting im- 
ages from the prior epoch(s). These components of the template 
image are astrometrically registered using catalog matching to 
align them with the WFC3-IR image of the current epoch. The 
astrometric registration for the SN search is done tile-by-tile 
and the output pixel grid is left in the natural unrotated frame 
of the observation. This contrasts with the CANDELS mosaic 
imaging pipeline (Koekemoer et al. 2011), which constructs a 
global astrometric solution across the whole field, and rotates 
every image to put north up and east to the left. These choices for 
the SN pipeline are designed to maximize the precision of the 
local inter-epoch registrations and to minimize dilution of the al- 
ready undersampled point-spread function (PSF) for single-visit 
drizzled images. 

Next, each template image is subtracted from the correspond- 
ing search epoch image, producing the difference images for SN 
discovery. Due to the very stable PSF of HST, the CANDELS 
images do not require any convolution with a PSF kernel to 
match conditions across epochs (Alard & Lupton 1998), as is 
commonly done in ground-based SN surveys. The CANDELS 
visits were constructed with small positioning shifts after each 
exposure, such that the two H band and two J band exposures 
together formed a four-point “box” dither pattern. This yields 
better sampling of the PSF and helps in the removal of detector 
artifacts from the final combined image. To take advantage of 
the full dither sequence, our SN searching was primarily done 
on a combined “/+//” image — simply the sum of the F125W 
and F160W difference images for each epoch. 

In the final stage of the data processing pipeline, we reprocess 
all the search epoch data, this time with fake SNe planted into 


the WFC3-IR survey images. These synthetic SNe enable a 
direct measurement of the detection efficiency of our human 
searchers (see Section 3). Each fake SN consists of a small 
image (~50 x 50 pixels) of a simulated point source, generated 
using the TinyTim software (Krist et al. 2011). The fake SN 
images are added to the WFC3-IR images at the FLT stage, 
after image calibration and before drizzling. These “faked” FLT 
files are then redrizzled, and the existing template images are 
subtracted off, resulting in a parallel set of “faked” difference 
images. 

2.2. SN Discovery 

To find SN candidates in the CANDELS WFC3-IR difference 
images, we used human searchers, who scanned each image 
by eye to detect significant deviations from the noise. We had 
~20 individuals regularly engaged in searching the CANDELS 
data, and searching tasks were assigned so that every WFC3-IR 
tile was examined by at least two people. Searchers recorded 
the position of all potential transient object detections, and 
assigned a quality grade. All transient sources that received 
a high- or moderate-quality grade were carefully vetted to pare 
down the list to transient sources that are very likely real SNe. 
The criteria for inclusion in this SN candidate list are: (1) a 
profile consistent with a point source, (2) detected in the J+H 
difference image and also individually in both the J and H bands, 
and (3) clean of evidence indicating that it could be a detector 
artifact (neighboring bad pixels, on the detector edge, etc.). For 
the rate analysis presented here, we also require that the object 
reached its peak magnitude in IR bands after the HST template 
images were collected (i.e., we reject any SN that has less flux 
in the search epoch than in the template epoch). 

Finally, we also discard a total of six objects that are positively 
classified as AGNs. These six were located at the center of a 
host galaxy that has observational indicators to classify it as an 
AGN (X-ray emission, spectral line broadening, prior optical/IR 
variability, etc.). Our final sample contains 65 SN candidates that 
meet these requirements. Table 3 lists the 14 SNe at redshifts 
z > L5, and Table 9 lists the remaining 51 SNe. In keeping 
with the practice of past HST SN surveys, we assign each SN 
a unique eight-digit name that indicates the field and the year 
of discovery, with the final three letters referencing our team’s 
internal “nickname” for each object. 28 Figure 2 shows “postage 
stamp images” with the detection images for the 14 SNe at 
z > 1.5, and the remainder are in the Appendix, in Figures 13 
and 14. 

2.3. Follow-up Observations 

Upon discovery, every SN was evaluated for possible follow- 
up observations with HST or ground-based telescopes. First, a 
redshift probability density function (PDF) was assigned, using 
pre-existing spectroscopy of the host galaxy when available 
and a photometric redshift (photo-z) when not. The photo-z 
estimates were derived from template fitting to the observed 
spectral energy distribution (SED) of each SN host galaxy 
(Dahlen et al. 2013). Then a preliminary SN classification (la or 
CC) was assigned by comparing the color and magnitude of the 
observed SN against a sample of synthetic SNe with redshifts 
drawn from the best available redshift PDF. These synthetic 
SNe were generated with the SuperNova ANAlysis software 
(SNANA; Kessler et al. 2009b) (see Section 4 for more details). 


27 


http ://www. stsci.edu/institute/ software_hardware/pyraf/stsdas 
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Presidents and other prominent figures from U.S. history. 
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GND13War,z= 1.689 


GSDllBus, z 


1.7 ± 0.1 


UPS 11 Hug, z= 1.76 ±0.03 



GND13Sto, z= 1.83 ±0.10 


GND12Bre, z= 1.88 


UDS10Wil,z= 1.914 




GND12Fai,z= 1.92 ±0.07 


GND12Col,z = 2.24 ±0.04 


EGSllTyl, z 


2.24 ±0.10 





GSD12Hum,z = 2.343 


I 



GSD12 


z= 2.37 


Figure 2. Detection images for 14 SN from the CANDELS fields with redshifts z > 1.5. Each image triplet shows //band (F160W) images with the template image 
on the left, the discovery epoch image in the middle, and the difference image on the right. All images have a width of about 6 arcsec, with north up and east to the 
left. The position of the SN is marked by (red) crosshairs in every frame. Discovery images for the other 51 SN with z < 1.5 are provided in Appendix B. 

(A color version of this figure is available in the online journal.) 


Any SNe with a redshift z > 1 and a color consistent with an 
SN la classification were then considered for possible follow- 
up with HST. Where necessary and whenever possible, the host 
galaxies of these high-priority targets were quickly observed 
(within ~1 week of discovery) with ToO spectroscopic obser- 
vations using ground-based observatories (primarily Gemini, 
Keck, and the Very Large Telescope (VLT)). The host galaxies of 
other SN candidates (CCSN and those with z < 1) were targeted 
for later spectroscopic observations from the ground to deter- 
mine precise redshifts, all reported in Table 4 (and in Table 10). 

Some of the most promising candidates for classification as 
SNe la at z > 1.5 were selected for supplementary imaging 
and/or grism spectroscopy with HST. Two of these, SNe 
GSDIOPri and UDSlOWil, have been presented elsewhere 
(Rodney et al. 2012; Jones et al. 2013). Due to the high cost 
of grism observations (at least 10 HST orbits are required to 
reach sufficient S/N in distant SNe), we applied strict criteria 
for selecting grism targets: (1) best available redshift z > 1 , 
preferably z > 1.5; (2) observed SN colors consistent with 
a (possibly reddened) SN la at that redshift; (3) observed SN 
magnitudes within ~1.5 mag of an SN la at that redshift (i.e., 
using a very weak prior around a standard ACDM cosmology); 
and (4) SN position allows for a grism observation without 
severe contamination. 

Without a slit to isolate the SN light in WFC3-IR grism 
spectroscopy, a high-z SN la candidate can most productively 


be observed if the trace of the SN spectrum can be positioned 
to avoid contamination from nearby galaxies. Thus, to satisfy 
the final criterion (4), the candidate must be well separated 
from the core of its host, or located in a host that is faint 
relative to the SN. We also require an orientation angle that 
avoids contamination of the SN spectral trace from the zeroth 
order and first order light of other nearby stars and galaxies. 
Of course, this orientation must also be accessible to HST at 
the time of observation, with suitable guide stars in range. In 
practice, these criteria were satisfied for only six SN candidates. 
The results of those observations are described in Section 6. 
Another 37 CANDELS SNe were followed with ToO imaging 
observations. These imaging targets included SN la candidates 
that satisfied some or all of the first three criteria, but were not 
suitable for grism observations, as well as some likely CCSNe 
that we were able to include in the same field of view as those 
primary targets. 

3. DETECTION EFFICIENCY 

Translating SN detections into an SN rate measurement 
requires characterization of the survey detection efficiency, i.e., 
the fraction of SNe that are detected by our human searchers. 
This recovery fraction is most strongly influenced by the S /N 
of the object in the WFC3-IR difference images. The SN host 
galaxy is also an important factor affecting SN detectability, as 
we discuss further in Section 3.1. 
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Table 3 

14 Supernovae with z > 1.5 (See Appendix B for the Remainder) 


Name 

R.A. 

Deck 

P(Ia|D,) a 

P (Ia|D h0St ) b 

ZSN C 

(±) 

z Source d 


(J2000) 

(J2000) 






COS12Car 

10:00:14.726 

+02:11:32.57 

0 62 +0 09 
u - oz -0.36 

0-80 

1.54 

(0.04) 

SN spec-z + SN phot-z 

GSDIOPri 

03:32:38.010 

-27:46:39.08 

1 00 +00 ° 
AUU -0.00 

1-00 + _ 0 0 % 

1.545 

(0.001) 

host+SN spec-z 

EGS13Rut 

14:20:48.106 

+53:04:22.12 

1 00 +00 ° 
iUU -0.00 

1 00 +00 ° 
A - uu -0.00 

1.614 

(0.005) 

host spec-z + SN phot-z 

GND13War 

12:36:54.761 

+62:12:16.70 

f) Q1 +0.02 
UU1 -0.01 

0 01 +00 ° 
UUi -0.00 

1.689 

(0.005) 

host spec-z 

GSDllBus 

03:32:42.776 

-27:48:07.10 

000 -°ri°oo 

0 00 +0 0 ° 
u - uu -0.00 

1.7 

(0.1) 

host+SN phot-z 

UDSllHug 

02:17:37.427 

-05:08:41.43 

0-82 !° 0 °2 5 i 

1 00 +00 ° 
A - uu -0.00 

1.761 

(0.025) 

host+SN phot-z 

GND13Sto 

12:37:16.778 

+62:16:41.43 

1 00 +00 ° 
A - uu -0.00 

1 00 +00 ° 
A - uu -0.00 

1.83 

(0.10) 

host+SN phot-z 

GND12Bre 

12:36:55.520 

+62:13:58.82 

0 00 +0 0 ° 
u - uu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

1.880 

(0.001) 

host spec-z 

UDSlOWil 

02:17:46.336 

-05:15:24.00 

1 00 +00 ° 
iUU -0.00 

1 00 +00 ° 
iUU -0.00 

1.914 

(0.001) 

host+SN spec-z 

GND12Fai 

12:36:15.822 

+62:15:56.50 

0 00 +0 0 ° 
u - uu -0.00 

0 00 +0 0 ° 
uuu -0.00 

1.92 

(0.07) 

host+SN phot-z 

GND12Col 

12:36:37.569 

+62:18:32.93 

1 00 +00 ° 
i>uu -0.01 

1 00 +00 ° 
iUU -0.00 

2.24 

(0.04) 

host+SN phot-z 

EGSllTyl 

14:20:12.944 

+52:57:10.60 

n 24 +013 

U - Z4 -0.15 

0 57 + 0 03 

U - J/ -0.04 

2.244 

(0.095) 

host+SN phot-z 

GSD12Hum 

03:32:15.500 

-27:50:50.02 

0 00 +0 0 ° 
u - uu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

2.343 

(0.001) 

host spec-z 

GSD12Qua 

03:32:11.723 

-27:49:11.72 

0 00 +0 0 ° 
uuu -0.00 

0 00 +0 0 ° 
uuu -0.00 

2.370 

(0.001) 

host spec-z 


Notes. 

a Type la SN classification probability from STARDUST, using the redshift-dependent class prior. Uncertainties reflect systematic biases due to the class prior and 
extinction assumptions (Sections 4.2 and 4.3). 

b Type la SN classification probability from STARDUST, using the galsnid host galaxy prior. Uncertainties reflect systematic biases due to the class prior and extinction 
assumptions. 

c Posterior redshift and uncertainty, as determined by the STARDUST light curve fit. 

d The host/SN values indicate whether the redshift is derived from the host galaxy, the SN itself, or a combination; spec-z/phot-z specify a spectroscopic or photometric 
redshift. A value of host+SN phot-z means the redshift is derived from a STARDUST light curve fit, with the host galaxy phot-z used as a prior. 


Table 4 

Host Galaxies of 14 Supernovae with z > 1.5 (See Appendix B for the Remainder) 


SN 

R.A. 

(J2000) 

Decl. 

(J2000) 

d 

(") 

d 

(kpc) a 

Morph b 

SED C 

£host 

(±) 

z Reference d 

COS12Car 










GSDIOPri 

03:32:37.991 

-27:46:38.69 

0.46 

9.6 

i 

SB 

1.545 

0.001 

Frederiksen et al. (2012) 

EGS13Rut 

14:20:48.113 

+53:04:22.07 

0.08 

1.7 

d 

A 

1.614 

0.005 

HST+WFC3 (A. Riess) 

GND13War 

12:36:54.787 

+62:12:16.60 

0.21 

4.3 

di 

SB 

1.689 

0.005 

HST+WFC3 (B. Weiner) 

GSDllBus 

03:32:42.776 

-27:48:07.10 

0.00 

0.0 

u 

A 

1.76 

0.53 

phot-z (T. Dahlen) 

UDSllHug 

02:17:37.415 

-05:08:41.53 

0.21 

4.2 

s 

P 

1.82 

0.13 

phot-z (T. Dahlen) 

GND13Sto 

12:37:16.823 

+62:16:42.65 

1.26 

25.5 

u 

A 

1.8 

1.2 

phot-z (T. Dahlen) 

GND12Bre 

12:36:55.520 

+62:13:58.79 

0.03 

0.6 

i 

SB 

1.880 

0.005 

Keck+MOSFIRE (J. Trump) 

UDSlOWil 

02:17:46.332 

-05:15:23.90 

0.12 

2.4 

s 

SB 

1.914 

0.001 

Jones et al. (2013) 

GND12Fai 

12:36:15.934 

+62:15:55.91 

0.98 

19.9 

sd 

SB 

1.77 

0.25 

phot-z (T. Dahlen) 

GND12Col 

12:36:37.514 

+62:18:32.66 

0.47 

9.6 

s 

A 

2.1 

0.2 

phot-z (T. Dahlen) 

EGSllTyl 

14:20:12.938 

+52:57:10.62 

0.06 

1.2 

sd 

SB 

1.95 

0.45 

phot-z (T. Dahlen) 

GSD12Hum 

03:32:15.585 

-27:50:50.43 

1.20 

24.3 

di 

SB 

2.343 

0.001 

Balestra et al. (2010) 

GSD12Qua 

03:32:11.713 

-27:49:11.29 

0.45 

9.3 

di 

SB 

2.370 

0.001 

VLT+Xshooter (J. Hjorth) 


Notes. 

a Physical separation between the SN and center of the host, computed from the measured angular separation in the preceding column, assuming a flat ACDM 
cosmology with Ho = 70, Q m =0.3. 

b Visual classifications for host galaxy morphology: s = spheroid, d = disk, i = irregular. 
c Template-matching classification of host galaxy SED: P = Passive, A = Active, SB = Starburst type. 

d Unpublished spectroscopic observations are given as Observatory+Instrument (j name of PI). Host galaxy photometric redshifts are marked as phot-z (T. Dahlen 2013, 
private communication). 


To measure our SN detection efficiency and explore the 
associated systematic biases, we generated a catalog of 2,000 
fake SNe. The catalog was drawn from a SNANA Monte Carlo 
simulation, such that the F160W magnitudes fill out a uniform 
distribution covering the range 21 < mu < 28, and the J — H 
colors were appropriate for Type la and Core Collapse SNe 
in the redshift range 0.1 < z < 2.8. Each fake SN was then 
assigned to a “host galaxy” drawn from catalogs of extended 


sources in the CANDELS fields. The separation from the host 
galaxy center for each fake SN was then selected randomly from 
a normal distribution centered on 0 with a standard deviation 
of 2 x R 5 q, where R 50 is the radius of an aperture containing 
50% of the host galaxy flux. This ensures that the fake SNe very 
roughly follow the distribution of host light (Kelly et al. 2008). 

With magnitudes, colors, and positions defined, we generated 
synthetic PSFs for each fake SN using TinyTim, and planted 
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Approx. Redshift for a SN la at Peak 
0.22 0.36 0.62 1.3 2.1 2.8 3.2 3.6 



Figure 3. SN detection efficiency measurements as a function of magnitude in 
the “J+H” band, taken as an average of the measured F125W and F160W 
magnitudes. Each point represents the fraction of fake SNe recovered by 
human searchers, with error bars indicating the standard deviation of the 
efficiency, computed using a Bayesian formalism (Patemo 2004). The best- 
fit model is shown as a solid (green) line, with best-fit parameters listed in the 
lower left. For reference, the equivalent best-fit curves for the J and H bands 
individually are shown in blue dotted and red dashed lines, respectively. The 
horizontal and vertical lines mark the 50% efficiency point for J+H detections: 
mso = 25.4 mag. The top axis marks the approximate redshift of a normal 
SN la with average extinction (Ay = 0.3) that would reach a peak brightness 
matching the J+H magnitude on the bottom axis. 

(A color version of this figure is available in the online journal.) 


them in the FLT images, as described in Section 2.1. As the 
searchers reviewed each difference image, they were unaware 
of the number, brightness, and location of the fake SN, so they 
recorded fake SN detections alongside detections of real SNe. 
After completing each search, the fake SN detections (and non- 
detections) were used to calculate the recovery fraction. 

Figure 3 shows the measured detection efficiency as a function 
of the “J+H' magnitude: the average of the F125W and 
F160W magnitudes. We fit the efficiency measurements with a 
functional form similar to that used by Sharon et al. (2007), but 
we use only a single parameter to characterize the exponential 
turnover, and we allow for the peak efficiency to plateau at a 
value less than unity: 

-l 

m — mso 

?7det(ra) = Y ] o x 1 + exp , (1) 

T 

where m is the apparent J+H magnitude, rj o is the maximum 
efficiency, mso is the magnitude at which the efficiency curve 
passes through the 50% line, and r characterizes the exponential 
roll-off. The best-fit curve shown in Figure 3 has mso = 25.4, 
r = 0.23, and rjo = 0.98. 

3.1. Missing SNe in Galaxy Cores 

One concern for systematic bias entering into these detection 
efficiency measurements is the possibility that many SNe are 
obscured by difference imaging artifacts in the cores of bright 



Figure 4. Fraction of CANDELS galaxies showing IR core residuals as a 
function of redshift, as determined from visual inspection of difference images 
generated by the SN data processing pipeline. 


galaxies. The shot noise in these bright pixels is naturally higher 
than in the outskirts, as photon counts are elevated in both 
the search epoch and the template. Additionally, minor cross- 
epoch registration errors can result in some residual flux in 
galaxy cores. In the CANDELS survey data these effects are 
both exacerbated by the under- sampled PSF of our single-epoch 
WFC3-IR images, as we have only two dithers per filter. 

As shown in Figure 3, we have measured our maximum 
detection efficiency r]o to be less than unity even for very bright 
SNe, due to the fake SNe that happen to land in the noisy cores 
of bright galaxies. Our SN rate measurements will therefore 
naturally account for a small fraction of SNe that are missed 
in this manner. However, this built-in correction is only valid if 
the distribution of positions for the fake SNe — relative to their 
host galaxy cores — is closely matched to the true distribution of 
the SN la population. Furthermore, it requires that the galaxies 
chosen for “hosting” our fake SNe are themselves representative 
of the population of SN la hosts. Our fake SN procedures were 
designed to meet these requirements at low and intermediate 
redshifts, but this does not necessarily carry over into the new 
high-z regime. 

To evaluate whether this effect might be introducing a 
strong bias at high-z, we visually inspected the CANDELS 
IR difference images and identified all galaxies that exhibited 
strong residuals. For each galaxy we tabulated the spectroscopic 
redshift or the best available photo-z from CANDELS catalogs. 
Comparing this redshift distribution for core residuals against 
the count of all galaxies as a function of redshift gives us a 
measure of the fraction of (detected) galaxies that might obscure 
SNe in their bright cores. As shown in Figure 4, the fraction is 
less than ~10% for all redshifts above 0.01, and less than ~2% 
for z > 1 — consistent with the value of rjo measured from fake 
SNe. This result suggests that any systematic bias from galaxy 
core residuals is very minor. Therefore, in the rate calculation 
we do not include any bias correction, and we do not add any 
contribution to the systematic uncertainty budget. 
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4. CLASSIFICATION 

To reach the final classification probabilities listed in Table 3 
(and Table 9), we used a Bayesian analysis of the observed multi- 
color light curves. This photometric classification approach 
was used for our full sample of 65 SNe, supplemented by 
spectroscopic evidence for six objects, as described in Section 6. 
An early version of this classifier was introduced in Jones et al. 
(2013) with the presentation of SN UDSlOWil. Here we will 
again briefly describe the classification procedure, emphasizing 
recent changes. 

4.1. The STARDUST Classifier 

Our photometric classification approach uses SNANA to 
generate simulations of SN la and CCSN light curves. The 
SN la simulations use the SALT2 model (Guy et al. 2010), 
which has free parameters for the date of peak (MJD p k), redshift 
(z), shape (jci), and color (C). The simulated CCSNe are drawn 
from the SNANA library of 42 CCSN light curve templates 
(26 Type II and 16 Type Ib/c). These templates are derived 
from the SN samples of the Sloan Digital Sky Survey (Frieman 
et al. 2008; Sako et al. 2008; D’ Andrea et al. 2010), Supernova 
Legacy Survey (Astier et al. 2006), and Carnegie Supernova 
Project (Hamuy et al. 2006; Stritzinger et al. 2009; Morrell 
2012). Each CCSN template defines the underlying shape and 
color of the synthetic light curves, which is then modified with 
free parameters for the date of peak, redshift, host extinction 
(Ay), and luminosity (Am, the shift in magnitudes relative to the 
peak of the assumed luminosity function). For this work, we fix 
the SALT2 model parameters a = 0.135 and = 4.1 (Scolnic 
et al. 2014), and for all simulated CCSNe we fix the extinction 
law to Ry = 3.1. 

Comparing these synthetic SNe to the observed light curves, 
we compute a likelihood using the x 2 statistic: P(D LC |0, la) oc 
exp(— x 2 /2), where the vector Dlc is the observed SN light 
curve and the vector 0 gives the parameter values for each re- 
alization of the SNANA models. We then apply priors for each 
model parameter (see Graur et al. (2014) for a detailed descrip- 
tion of these priors) as well as a redshift-dependent prior for 
the fraction of SNe that are Type la: P(Ia, z) (see Section 3 4.2 
below). Finally, we derive the total posterior probability that 
each object is an SN la, F(Ia|D L cX by marginalizing over the 
nuisance parameters, 0 , and applying Bayes’ theorem: 

T(Ia|D LC ) = P(\^z)P(Y) LC \0Aii)d0. (2) 

e 

The normalization factor k^c is defined by requiring that 
that posterior probabilities for all three primary SN classes 
(la, Ib/c, II) sum to unity. The custom-built software package 
that executes this procedure is named STARDUST: Supernova 
Taxonomy And Redshift Determination Using SNANA Tem- 
plates. The STARDUST code will be presented in full and pub- 
licly released in a subsequent paper (S. A. Rodney et al., in 
preparation). 

There are two notable differences between the STARDUST 
classification procedure applied here and that described in Jones 
et al. (2013). First, in this work we do not use a free parameter 
for flux scaling, 29 so the absolute values of the simulated 
SN fluxes are defined by the SN luminosity functions and 


29 The A in Equation (1) of Jones et al. 2013. 


cosmology (Q m = 0.3, Qa = 0.7, w = — 1 ) that are assumed 
in the SNANA simulations. To allow for some uncertainty in 
this baseline cosmology (or equivalently, introducing some 
increased scatter in the assumed SN luminosity functions) 
we include a non-zero model uncertainty term in the x 2 
calculation. 30 This is fixed at 8% of the simulated flux for SN la 
models and 10% for all CCSN models. Secondly, when the SN 
in question does not have a precise redshift from host galaxy 
spectroscopy, we use the host galaxy’s photometric redshift 
probability distribution (photo-z PDF) as the redshift prior. 

Column 4 of Table 3 (and Table 9) presents the final SN 
classification probabilities, which will be used in Section 7 for 
the SN la rate calculation. Figure 5 (and Figures 15-17) shows 
the maximum likelihood light curve fit for each SN, along with 
the associated best-fit model parameters. As described below, 
the systematic uncertainties associated with each classification 
probability are determined by varying two key priors that are 
not tightly constrained by observations: the assumed fraction 
of SNe that are of Type la and the distribution of host galaxy 
extinction. 

4.2. The Class Prior 

As with any Bayesian classification approach, the STAR- 
DUST classifier requires an input prior that quantifies the ex- 
pectation that any given SN is of Type la — before applying any 
information from the SN light curve. We first assume that our 
sample is composed entirely of “normal” SNe, meaning that 
we assume no contamination from any other transient sources. 
This is a fairly safe assumption: AGNs and variable stars are 
excluded by our discovery requirements, under- luminous SNe 
like the .la (Bildsten et al. 2007) or lax SNe (Foley et al. 2013) 
are well below our detection threshold, and super-luminous SNe 
(Gal- Yam 2012) have an intrinsic rate that is lower than that of 
normal SNe by a factor of about 10 4 (Quimby et al. 201 1). 

We then define a redshift-dependent class prior P(Ia, z) as the 
fraction of all normal SNe at any given redshift z that are Type 
la. Figure 6 shows the models used to define this prior and the 
associated systematic uncertainty. The baseline model (green 
curve) is anchored at z = 0 by the measured la fraction (Smartt 
et al. 2009; Li et al. 201 1), and then evolves at higher redshifts by 
following simple rate functions that match measured SN rates 
and theoretical expectations. 

A more complete statistical treatment would define a large 
number of plausible models for the fraction of SNe that are 
SNe la, assigning each an appropriate weight based on current 
observations, and then marginalize over those many discrete 
priors to get a posterior probability that is not uniquely guided 
by the single choice of a baseline model. That approach is 
computationally expensive and will require further refinement 
of the STARDUST classifier. For this work, we have chosen 
to treat this choice of prior as a component of our systematic 
uncertainty budget. We take the baseline prior described above 
as our mid-rate model and then define two more models, labeled 
the high- and low-rate priors. These two respectively maximize 
and minimize the fraction of SNe that are assumed to be of 
Type la at any given redshift, and are shown in Figure 6. These 
bounding models offer a conservative estimate of the systematic 
uncertainty, because they are at the extreme limit of plausibility 
(if either were correct it would imply that the constraints from 
past rate measurements were all systematically wrong by more 
than 2a). 


30 a? in Jones et al. 2013. 
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COS12Car 
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Figure 5. Light curves of 14 SN with z > 1.5. Each panel shows the observed photometry (in Vega mags) from CANDELS imaging with filled points and arrows 
for 3 a upper limits. Error bars are typically less than the size of the points. Curves depict the maximum probability light curve fit for the most probable SN class, as 
determined by STARDUST. Classification probabilities and redshift as determined by STARDUST are listed on the right side of each panel, along with the parameters 
of the single model that delivers the highest posterior probability. Light curves for SN at z < 1.5 are shown in Appendix B. 


One might be concerned about the apparent circularity of 
using a redshift-dependent P( la, z) prior based on measured 
SN la rates in the service of a new SN la rate measurement. 
However, the bounding assumptions for our classification prior 
should ensure that our systematic uncertainty estimates account 
for this. To test that assertion, in Appendix A we evaluate 
an alternative prior that is based on the SN host galaxies 
and does not evolve with redshift. Tables 3 and 9 record the 
resulting STARDUST classifications using this modified prior 
as P (la | D host ) in Column 5. Table 8 reports the final effect of this 
prior switch on the observed count of SNe la and the volumetric 
rates. 


4.3. Host Ay Distribution 

Another prior that can strongly affect the final classification 
probabilities is the assumed distribution of host galaxy extinc- 
tions, P (Ay). As with the SN la fraction prior, we employ a 
baseline assumption (our mid-dust model) and two bounding 
assumptions (high-dust and low-dust) to constrain the possible 
systematic bias. 

In keeping with observations, our dust models assume that 
the CCSN population suffers from significantly more dust 
extinction than the SN la population, at all redshifts (e.g., Smartt 
et al. 2009; Drout et al. 2011; Kiewe et al. 2012; Mattila 
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Figure 6. Deriving the redshift-dependent prior for SN class fractions. The 
top panel shows SN la rates and the middle panel shows CCSN rates. Both 
have observed rates plotted as open symbols. In the la case these are average 
values from all non-redundant field SN surveys. The CCSN rate points are the 
collection from Dahlen et al. (2012). In each panel the overlaid solid lines show 
three versions of a simple empirical model for the SN rates, and in the CC 
panel the magenta line traces the cosmic star formation history (see the text for 
details). The bottom panel plots the fraction of all SN explosions that are of Type 
la, derived from pairs of curves drawn from the top two panels and anchored 
to 25 ± 5% at z = 0 (Smartt et al. 2009; Li et al. 2011). These relative rate 
assumptions provide the high-, mid-, and low-rate priors that are used to derive 
classification probabilities and associated uncertainties for observed SNe. 

(A color version of this figure is available in the online journal.) 

et al. 2012). Our three dust models are generated from the 
positive half of a Gaussian distribution centered at Ay = 0 
with dispersion cr , plus an exponential distribution of the form 
e ~ Ay/ 1 . The parameter Ro gives the ratio of the height of the 
Gaussian to the height of the exponential, at Ay = 0. The 
defining parameters and the expectation values for these three 
distributions are summarized in Figure 7. 

When simulating SN la with the S ALT2 model, the SN color 
is defined by the SALT2 C parameter. This color term comprises 
both the intrinsic SN color as well as reddening from host galaxy 
dust. The distribution of C values can therefore be described as a 
convolution between a narrow Gaussian (the intrinsic dispersion 
of SN la colors) and a function describing the distribution of 
host galaxy extinctions. Following Barbary et al. (2012) and 
Scolnic et al. (2014), we approximate the Ay distributions of 
previous SN la studies by modifying the red side of the SALT2 
C distribution so that the simulated SN colors match the output 
of that convolution. 

Specifically, our high-dust model for SN la matches the 
baseline Ay distribution used by Neill et al. (2006): a Gaussian 
with cr = 0.62. The mid-dust model is equivalent to the 
exponential distribution of Kessler et al. (2009a): P(A V ) = 
^-Ay/r, w j t | 1 r _ o.33. Our low-dust model for SN la assumes 
minimal dust extinction, using a narrow Gaussian with a = 0.15 
plus a shallow exponential with r = 0.15. A more complete 


SALT2 C 
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host extinction ( Ay ) 

Figure 7. Prior probability distributions for the SN host galaxy extinction, as 
used in the STARDUST classifier code. The top panel shows the three priors 
applied to the SN la models, and the bottom panel shows the equivalent priors 
used for CCSN models. In both cases the high-dust model is shown as a (red) 
dashed line, the mid-dust model as a (green) solid line, and the low-dust model 
as a (blue) dash-dotted line. Each model is composed of a Gaussian and/or 
an exponential function (see the text for details) and the parameters for those 
components are listed in the legends. Each curve is labeled with its expectation 
value (Ay), giving the “weighted average” of the host galaxy extinction for that 
model. 

(A color version of this figure is available in the online journal.) 

treatment of host galaxy dust would include a prescription for 
the redshift dependence of these extinction distributions, as 
SN hosts are expected to be dustier at redshifts approaching 
z = 2 (Mannucci et al. 2007). As we will see in Section 7, this 
systematic uncertainty is not a dominant component of the error 
budget, so redshift dependence is left for future work. 

To determine the combined systematic effects from the 
SN la fraction prior and the dust assumptions, we compute 
each SN classification probability nine times: three for each 
SN rate prior x three for each dust model. The mid-rate + mid- 
dust combination gives us our baseline classification probability, 
which dictates how much each individual SN contributes to the 
total count of observed SN la, A 0 b S . The extrema from this set 
of nine probabilities then provide the systematic classification 
errors, which propagate directly into the systematic uncertainty 
on the SN la rate. 

4.4. STARDUST Validation Test 

A full investigation of the accuracy of the STARDUST 
classification code is beyond the scope of this paper. It is useful, 
however, to examine a simple validation test to demonstrate that 
this classifier is not grossly biased or ineffective. To that end, 
we have applied the STARDUST classifier to the “Gold” sample 
of 31 SNe from the GOODS and PANS surveys (Strolger et al. 
2004; Riess et al. 2007) that have spectroscopic classifications. 
These surveys were carried out using the HST Advanced Camera 
for Surveys (ACS), and share many of the survey design 
characteristics of the CANDELS SN program. STARDUST 
correctly classifies 29 of the 31 SNe (93.5%), using only their 
redshifts and photometric data. This demonstrates that we have 
a low false negative rate with STARDUST, i.e., we rarely 
misclassify a true SN la. Unfortunately, this validation test is 
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not sensitive to false positives — true CCSNe being misclassified 
as Type la — because we only have a single spectroscopically 
confirmed CCSN in this Gold sample. Preliminary testing of 
the STARDUST classifier using simulated SNe suggests that 
the la sample purity for photometrically classified SNe could be 
on the order of 95% (these validation tests will be presented in 
a future paper). 

5. HOST GALAXIES 

Host galaxy information is recorded in Table 4 for the SNe at 
z > 1.5, and in Table 10 for the low redshift SNe. As can be seen 
in Figure 2 (and Figures 13-14), most of the 65 CANDELS SNe 
can be unambiguously associated with a host galaxy, because the 
host is isolated, or the SN is clearly embedded within the stellar 
light of a single galaxy. There are, however, a few exceptions. 

For SN COS12Her, there are two host galaxy candidates: 
the nearest and brightest has a photometric redshift of 0.403 
to 0 /! , but the observed SN colors cannot be adequately matched 
by any normal SN template in that redshift range. The second 
COS12Her host candidate has a photo-z = 1 . 1 0 +q \ 6 9 . At this 
higher redshift, the STARDUST classifier finds a very good 
match to the observed light curve with a Type II-P template. 

SN GND13Sto is separated by several arcseconds from all 
nearby galaxies. Of the six galaxies within 5 arcsec of the SN 
position, five have a photo-z distribution that peaks close to 
z = 1.8, including one with a spectroscopic redshift from the 
Spitzer Infrared Spectrograph of z = 1 .80 zb 0.02 (Murphy et al. 
2009). This is suggestive of a small cluster or group of galaxies 
at that redshift, with SN GND13Sto possibly associated with a 
low surface brightness group member or tidal stream. Indeed, 
applying STARDUST to the well-sampled SN GND13Sto light 
curve (and allowing for a broad redshift range, z = 1.8 zb 1.2), 
we find the maximum likelihood match is an SN la template at 
z = 1.86 =b 0.05. 

There are four SNe for which the host galaxy is barely de- 
tectable in the deep IR imaging mosaics from CANDELS. These 
are SNe GND12Kin, EGSllNix, GND13Gar, and GSDllBus. 
All four of these objects lack a clear spectroscopic redshift from 
their host, so we are limited to using photometric redshifts for the 
STARDUST priors. In all of these cases, with STARDUST we 
find good template matches within the allowed redshift range. 

5.7. Morphology and SED Type 

The SN host galaxies (along with all CANDELS galaxies) 
were classified visually by members of the CANDELS team into 
three morphological categories: spheroid, disk, and irregular. 
Visual classifications were performed using template images so 
that the presence of the SN did not bias the classification. Each 
galaxy can be assigned to multiple categories, so we also include 
two intermediate categories: spheroid+disk and disk+irregular. 
These morphological classes roughly correspond to broad bins 
over the Hubble sequence. This is appropriate for classifying 
galaxies at high redshift where distinguishing between, say, an 
E and an SO galaxy is more difficult and less meaningful. Details 
of the CANDELS morphological classification procedure are 
presented in Kartaltepe et al. (2012). 

We also record the “SED type” for each SN host galaxy, 
determined by matching the full galaxy SED against a set of 
templates, using the GOODZ code (Dahlen et al. 2010). The 
GOODZ template library is segregated into three groups, labeled 
according to the amount of ongoing star formation: passive 
(early type), active (late type), and starburst. We use the best- 


matching SED template for each CANDELS SN host galaxy to 
assign it to one of those bins. 

For two of the SNe (COS12Car and GND12Daw) there is 
no discernible host at the location of the SN and no nearby 
galaxy presents a plausible host candidate. For both of these 
objects we do have spectroscopic redshift information from 
the SNe themselves, as detailed in Section 6. For the other 63 
SNe in our sample, 10 host galaxies are classified as spheroids, 
15 as spheroid+disk, 17 as disk, 7 as disk+irregular, and 8 as 
irregular. For six of our SNe, the host galaxy is detected, but 
is too faint for reliable visual classification, so we report the 
host morphology as “unclassifiable.” For the 63 objects with 
detectable host galaxies we have 2 passive, 24 active, and 37 
starburst-like SEDs. 

6. GRISM SPECTROSCOPY 

There are six objects in our sample for which we collected 
useful HST grism spectroscopy of the SNe themselves. SN 
GSDIOPri, a Type la SN at z = 1.55, was presented in Rodney 
et al. (2012) with an analysis of the host galaxy in Frederik- 
sen et al. (2012). SN UDSllWil, a Type la SN at z = 1.91, 
was described in Jones et al. (2013). Figure 8 presents grism 
spectra for the remaining four: SNe GSDllWas, GND12Daw, 
GND13Gar, and COS12Car. In all of these cases the signal-to- 
noise ratios and rest-frame wavelength coverages are insuffi- 
cient for a purely spectroscopic classification. Rather, as with 
GSDIOPri and UDSllWil, we used the spectroscopic infor- 
mation to supplement the STARDUST photometric classifier, 
leading to a more robust classification. 

The host galaxy of SN GSD1 lWas has a photometric redshift 
of z = 1.04 zb 0.3. We obtained a spectrum of SN GSDllWas 
with the WFC3-G141 grism, shown in Figure 8 (top left). Here 
we can see hints of an absorption feature at ^14000 A. At 
a redshift of z ~ 1.3 this feature can be explained as the 
characteristic Si n absorption trough seen at rest-frame ~6150 A 
in SN la spectra. Photometric classification of this SN with 
STARDUST agrees, finding the object is best matched by an 
SN la template at z = 1.3 zb 0.05. (see the light curve plot in 
Appendix B, Figure 17). 

SNe GND12Daw, GND12Gar, and COS12Car all have no 
detectable host galaxy in any optical or NIR band, and no 
neighboring galaxies have redshifts that allow for acceptable 
light curve template matches in STARDUST. The most likely 
explanation is that these SNe reside in very low surface 
brightness galaxies, too faint for detection even in our deep 
HST imaging. 

The spectrum for GND12Daw shows hints of a broad emis- 
sion feature at 12000 A, which could be Ha emission, if the ob- 
ject is at z = 0.830. This could be interpreted as strong Balmer 
line emission from an otherwise very faint host galaxy, or it 
could be showing the Ha emission from the SN itself — charac- 
teristic of Type II-P spectra. Given the very low signal-to-noise 
ratio in this spectrum (it was derived from just a single orbit 
of HST observations), this alone would be weak support. How- 
ever, when allowing STARDUST to search over a redshift range 
0.1 < z < 2.0, we find that a Type II-P light curve template con- 
sistently provides the strongest match to the broad light curve 
shape of this SN, and the solution at z ~ 0.8 provides > 90% 
of the total likelihood. 

For GND12Gar (upper right) the absorption trough at 
~7700 A provides a key observable that can anchor the fit 
and define the age of the SN. If this feature corresponds to 
Can absorption, then that would fix the object’s redshift to 
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Figure 8. Observed spectra and spectral template matches for four CANDELS SNe. The spectra for GSD1 lWas (top left) and GND13Daw (top right) were collected 
using the G141 grism on HSTs WFC3-IR detector. The GND13Gar spectrum (lower left) used the ACS G800L grism, and the spectrum of COS12Car (lower right) 
combines observations from both the G102 and G141 WFC3 grisms. In each panel the observed flux is shown in gray, with binned points overlaid in blue, and the 
best-fitting template spectrum in red. 


z — 1.07 =t 0.02. At this redshift the light curve is matched very 
well by SN la templates, and no other redshift or SN class can 
provide a better light curve match. 

The strongest spectral constraint for SN classification comes 
from SN COS12Car. For this object we have observations with 
both the WFC3 G102 and G141 grisms. Fitting the composite 
spectrum with the SuperNova IDentification (SNID) software 
(Blondin & Tonry 2007), we find the best template match is a 
Type la SN at z = 1.59. Once again we find that the STARDUST 
photometric classification agrees well with this spectroscopic 
information: an SN la light curve at z ~ 1.6 provides the best 
available light curve template. 


7. THE VOLUMETRIC SN la RATE 

To convert the observed SN counts into a volumetric rate, 
we use an approach similar to Dahlen et al. (2008) and Rodney 
& Tonry (2010). We first divide up the detected SNe into four 
redshift bins of width A z = 0.5. The total contribution to the 
SN la count from each observed SN is equal to the la classi- 
fication probability for that object. For objects with uncertain 
redshifts, this fractional contribution is distributed over multi- 
ple redshift bins according to the integrated area of the redshift 
PDF. Adding up all the fractional counts gives us the total ob- 
served SN la count as a function of redshift: A 0 b s (z). Statistical 
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Table 5 

Observed SN la Counts and the Volumetric Rate 


Redshift 


Observed Count a 


Control Count b 


SN Rate c 



Vobs 

^ VAiss 

^-Vsyst 

Vtrl 

<5 N syst 

SNR 

<$SNR sta t 

<5SNR syst 

0.25 

1.46 

+ 2.46 

- 1.06 

+ 0.48 

- 1.44 

4.10 

+ 0.01 

- 0.04 

0.36 

+ 0.60 

- 0.26 

+ 0.12 

- 0.35 

0.75 

7.19 

+ 3.80 

- 2.62 

+ 1.94 

- 2.54 

14.11 

+ 0.53 

- 1.71 

0.51 

+ 0.27 

- 0.19 

+ 0.23 

- 0.19 

1.25 

8.47 

+ 4.02 

- 2.85 

+ 0.45 

- 2.04 

13.16 

+ 2.46 

- 4.05 

0.64 

+ 0.31 

- 0.22 

+ 0.34 

- 0.23 

1.75 

5.54 

+ 3.49 

- 2.28 

+ 0.17 

- 0.61 

7.67 

+ 3.47 

- 2.99 

0.72 

+ 0.45 

- 0.30 

+ 0.50 

- 0.28 

2.25 

1.24 

+ 2.39 

- 0.96 

+ 0.13 

- 0.16 

2.52 

+ 1.72 

- 1.07 

0.49 

+ 0.95 

- 0.38 

+ 0.45 

- 0.24 


Notes. 

a Statistical uncertainties reflect the limits that contain 68% of the Poisson distribution. Systematic uncertainties are due to the 
assumed dust model and rates prior. 

b Systematic uncertainties are due to the assumed dust model. 
c The SN la rate measurements in units of SNuVol = 10 -4 yr -1 Mpc -3 /i 7 o 3 . 


uncertainties for each bin are defined by the points encompass- 
ing the central 68% of the Poisson distribution. 

We then use Monte Carlo simulations to compute a “con- 
trol count” for each bin, N ctr i(z), which is the expected num- 
ber of SNe la that would be detected if the cosmic SN la 
explosion rate were constant for all redshifts at 1 SNuVol = 
10 -4 yr -1 Mpc -3 h^ 3 . By simulating SN la light curves within 
the context of the CANDELS survey, the computation of N ctr \(z) 
incorporates both the survey volume and the control time (the 
time interval over which any given SN is visible to our survey). 

We again use SNANA as our simulation engine, this time 
generating 100,000 SNe la based on the SALT2 light curve 
model (Guy et al. 2010). The light curve for each synthetic SN 
is determined by a set of four variables: date of peak brightness 
mjd pk, redshift z, SALT2 shape parameter x\, and SALT2 color 
parameter C. The C parameter in S ALT2 includes both intrinsic 
SN color as well as reddening from host galaxy dust. Each 
redshift z is drawn from the range 0 < z < 2.5, following the 
constant volumetric rate assumption. To translate this redshift 
into a luminosity distance, we use our baseline cosmology: 
Q m = 0.3, Qa = 0.7, w = — 1, Hq = 70. Values for x\ 
are drawn from a normal distribution with mean and dispersion 
from Kessler et al. (2009a): x\ = —0.13, a Xl = 1.24. The color 
parameters are draw from a bifurcated Gaussian distribution 
with C = 0.04, <7 q = 0.08, and ctq = 0.25 — parameters that 
match the mid-dust model described in Section 4.3 and Figure 7. 

To choose values for mjd v k we first establish the width of the 
survey window at any given redshift. By examining simulated 
SN la light curves in F125W and F160W, we find f m i n and f max , 
the minimum and maximum dates relative to peak for which 
each SN would be detectable to our survey. Here, detectability 
is defined by measuring the change in flux relative to a template 
epoch 52 days prior, and requiring that the corresponding J+H 
magnitude is brighter than the 50% detection threshold seen in 
Figure 3. The allowed range for the simulated mjd pk values at 
redshift z is then [MJD first - Ufe), MJDi ast + t min (z)], where 
MJDfi rst and MJDi ast are the epochs for the first and last search 
epoch, respectively. For each redshift, random mjd v k values 
are then drawn from a flat distribution spanning this survey 
window. 

Each synthetic SN is “observed” in the SNANA simulator 
using survey parameters that match the actual operations of the 
CANDEFS program, as given in Tables 1 and 2. For the wide 
fields (COSMOS, EGS, UDS, and the wings of the GOODS 
fields) we only have a single search epoch, so these fields are 
simulated together as the “CANDELS -Wide” search field. The 


10 epoch GOODS-S and GOODS-N Deep fields are treated 
separately, but all observational parameters are computed in the 
same way. Due to the excellent stability of the HST photometric 
system, we adopt a single set of average values for zero points 
and detector noise. The total area in each field reflects the area 
in which SN searching can be done, i.e., the area covered by 
the SN search epoch and at least one prior epoch. The cadence 
between epochs is nominally 52 days, but the actual separation 
in time varies from pointing to pointing due to HST scheduling 
constraints. For this simulation we use a mean cadence for each 
field and each epoch, weighted by the area available for SN 
searching. Finally, we use the detection efficiency curve of 
Figure 3 to define the probability of “detecting” each simulated 
SN in any given epoch. 

Counting the number of detected synthetic SNe in each 
redshift bin gives us the control count, which carries units of 
SNuVol -1 . The observed volumetric rate of SN la explosions is 
simply the ratio 


SNR Ia (z) = 


Nobs (z) 
N c trl(z) 


( 3 ) 


The measured values for N 0 b s (z), N ctr i(z), and SNR(z) from 
the CANDELS survey are given in Table 5 along with uncer- 
tainty estimates due to statistical noise (Poisson errors) and 
systematic biases. The total sample size is quite small, with 
only ~21 SNe la across all redshifts and fewer than seven in 
each bin. This means that the statistical errors are substantial, 
roughly equal to, or greater than the systematic uncertainties in 
every redshift bin. One cannot infer a clear trend with redshift 
from these data alone, but rather we must evaluate them within 
the context of other rates measurements and SN la progenitor 
models. 


7.7. Systematic Uncertainties 

In preceding sections we considered three principal sources of 
systematic biases: (1) missing SN detections due to subtraction 
artifacts in the cores of bright galaxies, (2) the assumed fraction 
of SNe that are of Type la as a function of redshift, and (3) 
the assumed distribution of host galaxy dust extinction values. 
We have determined that bias from the first source is negligible. 
The second is examined in more detail in Appendix A, and is 
reflected in the systematic uncertainty estimates for the count 
of observed SNe la (Column 4 of Table 5). The third item also 
affects the control count (Column 6). 

Other potential sources of systematic bias include: errors in 
the luminosity functions for SN subclasses, biases in the SN la 
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Figure 9. Measured SN la rates from CANDELS and other surveys. The CANDELS volumetric SN la rate measurements are shown as large orange squares, spanning 
five redshift bins of width dz = 0.5. For these CANDELS points, the systematic uncertainties are shown as broad orange bars, while the thin vertical error bars 
show the combined systematic and statistical uncertainty. Four other high redshift SN surveys are highlighted: gray circles for the Subaru Deep Field (SDF; Graur 
et al. 201 1), blue downward triangles for volumetric (not cluster) rates from the Cluster Supernova Survey (CSS; Barbary et al. 2012), green upward triangles for the 
GOODS and PANS surveys (Dahlen et al. 2008), and magenta diamonds for the CLASH survey (Graur et al. 2014). Assorted ground-based surveys are plotted as gray 
circles, as in Figure 1. 


model or the CCSN template libraries, and peculiar detection 
biases from individual human searchers. For this work, these 
contributions to the systematic error budget are assumed to be 
insignificant. Future analysis with the full CANDELS+CLASH 
sample will revisit this assumption and evaluate these systematic 
error sources. 

8. TESTING SN la PROGENITOR MODELS 

The measurement of volumetric SN la rates at high redshift is 
principally motivated by two astrophysical investigations. First, 
it directly informs our understanding of the cosmic enrichment 
history, as SNe la are a primary source for Fe group elements in 
the universe (e.g., Wiersma et al. 2011). Second, by measuring 
the delay between star formation and SN explosion through the 
DTD formalism, one can draw inferences about the nature of 
SN la progenitor systems. In this work we limit our discussion 
to the latter, beginning with a comparison to other published 
SN la rates, then evaluating new constraints on progenitor mod- 
els, and finally making a projection toward future improvements. 

8.1. Comparison to Earlier Rate Measurements 

Figure 9 presents the CANDELS SN la rates within the 
context of other rate measurements from the literature. The 
CANDELS rate measurements are shown in five bins of width 
A z = 0.5 as large orange squares. Rate measurements from 13 
ground-based surveys are plotted as small gray circles, reaching 
out to z = 1.1. Four surveys that previously extended the rate 
measurements to z ~ 1.5 are highlighted with larger colored 
points (see the caption for details). 

As with past HST surveys, our survey volume is too small 
to add any useful new information at z < 1 , but the general 


agreement with ground-based surveys is an important validation 
that our rate measurements are realistic. For a more informative 
comparison, we turn now to the high-z side of the plot. 

Before CANDELS and CLASH, there were just three surveys 
with any SN la rate measurements above z~ 1.2. First were 
the GOODS SN surveys, which used the HST ACS camera 
to measure the SN la rate to z ~ 1.8 (GOODS; Strolger 
et al. 2004; Dahlen et al. 2004). These data were interpreted 
as showing a flattening or a downturn in the SNR(z) at z > 1.2, 
a trend that garnered support from additional HST observations 
and independent analyses (Dahlen et al. 2008; Kuznetsova 
et al. 2008), including another HST+ACS program, the Cluster 
Supernova Survey (CSS; Barbary et al. 2012). 31 Ground-based 
rate measurements from the Subaru Deep Field (SDF) survey 
also reached out to z ~ 1.8, though these were more susceptible 
to systematic biases due to the absence of light curve information 
and spectroscopy for SN classification (Poznanski et al. 2007; 
Graur et al. 2011). As the CANDELS and CLASH surveys 
began, it was still an open question as to whether we had now 
seen the peak in the SN la rate, or if it was continuing to rise 
beyond z ~ 1.2. 

The GOODS, CSS, and SDF surveys all used optical bands 
that correspond to rest-frame near-UV wavelengths at high 
redshift. For an SN at redshift z ~ 1.5, observations in the z 
band (~9000 A) are sampling the rest-frame U band, while the 
observer’s i band (~8000 A) reaches well into the rest-frame 
near-UV. At these wavelengths the available SN light curve 
templates for use in photometric classification are poor, because 
most nearby SN surveys do not observe SNe in the UV. Both 


31 This was a survey of galaxy clusters, but the work of Barbary et al. (2012) 
presented volumetric SN la rates from the SNe detected outside the clusters. 
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SNe la and CCSNe also exhibit more natural heterogeneity 
at these blue wavelengths, and this is all compounded by 
a greater sensitivity to dust obscuration in the UV. Thus, 
optical-wavelength surveys were more susceptible to both of 
the components that dominate the systematic uncertainties of 
high-z SN la rate measurements: classification bias and dust 
obscuration. By contrast, the CANDELS survey utilized J and H 
infrared bands that sample rest-frame optical wavelengths, even 
out to redshift z ~ 2.5. The CANDELS rates should therefore 
be less strongly affected by those systematic biases. 

At z — 1.25 the CANDELS rate is substantially lower than 
all past measurements, though still consistent at the l-2cr level. 
The CANDELS rate then climbs slightly in the bin at z = 1.75, 
where it is completely consistent with past measurements. 
CANDELS is the only survey with any detections beyond z = 2, 
and there we have only a single object with a strong probability 
of being an SN la (SN GND12Col in the GOODS-N field, at 
z ~ 2.24). The rate formally shows a decline to z = 2.25, 
although this change is much smaller than the uncertainties. 
The CANDELS rates are fully consistent at all redshifts with 
the similarly derived rates from CLASH, which are also quite 
low relative to past surveys (Graur et al. 2014). 

Due to the small sample sizes and large uncertainties, none 
of these individual high-z surveys has sufficient precision to 
clearly delineate the shape of the SNR(z) curve. From Figure 9 
we can only say that the SN la rate rises steadily to z ~ 1.2, and 
then is flat or slowly declining at redshifts z > 1.2. 

Each independent analysis of SN la rates makes slightly 
different assumptions about host galaxy extinction and each 
takes a different approach to SN classification. These differences 
become particularly important at z > 1 where observed rates are 
dominated by HST SN surveys, which have much less complete 
spectroscopic information. Here the potential for systematic 
biases is much greater, as a larger fraction of SN classifications 
and host galaxy redshifts rely on photometric data alone. 

An optimal approach would be to treat the past decade of HST 
SN surveys as a single composite survey. One could compute 
the rates from all the HST SN surveys together, using the same 
SN classifier(s), consistent models for (redshift-dependent) host 
galaxy extinction, and the best available host galaxy redshift 
information. Such an effort is beyond the scope of this work, 
but will be an important contribution for future DTD tests. 


Age of Universe [Gyr] 

13.0 3.0 1.0 0.5 



Figure 10. Cosmic star formation rate (CSFR) as a function of redshift. Points 
show the compilation of recent CSFR measurements from Behroozi et al. 
(2013), adopting from those authors the corrections for dust attenuation and 
more realistic systematic errors. The solid line shows the best- fit double power 
law model from Behroozi et al. (2013), and the shaded region demarcates the 
la systematic uncertainties. 


Belczynski et al. 2005): 


0 

SNR(f) = ^Kr lu ^ fp 


for t < 0.04 Gyr, 
for 0.04 < t < 0.5 Gyr, 
for t > 0.5 Gyr. 


( 4 ) 


Here ^i a indicates the efficiency of generating SN la progenitor 
systems, in units of SN la, yr -1 M 0 _1 , and f P sets the fraction 
of all SNe la that arise from the prompt channel. The constant 
K is defined by the time thresholds that delineate this model: 


K — lu(^max/0 )/(0 ^min)> (5) 


8.2. Isolating the Prompt SN la Fraction 

To examine how the observed SN la rate can inform the 
modeling of SN la progenitors, we will employ a simple 
toy model, motivated by a variety of recent observations and 
theoretical predictions. For a complementary analysis using 
DTD predictions from binary population synthesis modeling, 
see Graur et al. (2014). Multiple lines of evidence now suggest 
that the overall shape of the SN la DTD follows a t~ l power 
law for times t > 500 Myr (see Maoz & Mannucci 2012, for a 
recent review). At short delays, t < 500 Myr, the evidence is 
much less definitive, and this is the region where the CANDELS 
observations may provide unique new insight. Thus, our primary 
question is: what fraction of SNe la explode within 500 Myr of 
their formation? 

To isolate this “prompt SN la fraction,” we define a bifurcated 
DTD model: the long-delay component follows a t~ l distribu- 
tion for all times t > 500 Myr, and the prompt component is set 
to be constant with time for t < 500 Myr, down to a lower limit 
of r min = 40 Myr (the shortest possible time to reach explosion; 


where t m i n = 0.04 Gyr as defined above, t\ =0.5 Gyr marks the 
abrupt transition from the constant rate to the power law, and 
t max = 13.3 Gyr is the maximum age of a WD in the current 
universe — using our assumed ACDM cosmology and assuming 
star formation began at z = 20. For these values, we have 
K = 7. 132. With this simple DTD model, we can allow rji a and 
fp to be free parameters, and fit to the data to find the observed 
efficiency and prompt la fraction. 

To convert from this DTD model into a prediction for SN la 
rates, we convolve this DTD with a parameterized representation 
of the cosmic star formation history, giving us a prediction 
for the observable SNR(z). For this exercise we use the recent 
compilation of measurements of the cosmic star formation rate 
(CSFR(z)) from Behroozi et al. (2013), shown in Figure 10. 
The precise shape of the CSFR curve at z > 2 is still 
a matter of debate, but for our purposes here we take the 
Behroozi et al. curve and associated systematic uncertainties 
to be representative of the current state of the art (but see 
Graur et al. 2014, for further evaluation of star formation history 
variation). 


15 


The Astronomical Journal, 148:13 (28pp), 2014 July 


Rodney et al. 


Table 6 

Observational Constraints on the SN la Delay Time Distribution 


Sample 

*7la 

fp 

M a/M* b ’ c 

N Ia/MA d 


(10 -4 Mq~ 1 yr -1 ) 


( 1 0 — 3 Mq -1 ) 

(10- 3 Mq- 1 ) 

CANDELS & CLASH 

9 9 c +1.36 +0.72 
-1.18 -0.15 

rj 91 +0.34 +0.49 
U - Zi - 0.21 - 0.12 

n 70 +1.52 +2.11 
U - /y -0.50 -0.16 

0 6n +0.97 +0.87 
uou -0.37 -0.46 

Ground 

1 78 +0.24 +0.43 
-0.23 -0.21 

0 SO +0.09 +0.05 
-0.10 -0.04 

O 8 zL +0.39 +0.44 
11.84 _ Q 2? _ 0 54 


All 

1 zrrj +0.24 +0.25 
iDU -0.23 -0.59 

n ^7 +0.09 +o.io 
- 0.10 -0.26 

O QR +0.43 +0.87 
-0.30 -0.46 

n 7Q +0.88 +1.09 
u -0.72 -0.75 


Notes. 

a Errors give first statistical then systematic uncertainties. 
b Assuming a t~ l delay time model of the form given in Equation (4). 

c Using the Behroozi et al. (2013) cosmic star formation history, which assumes a Chabrier (2003) IMF. 
d Using the SN la rate data directly, without any DTD model assumption. 


The construction of our bifurcated DTD model is reminiscent 
of the two-component “A+B” model (Mannucci et al. 2005; 
Scannapieco & Bildsten 2005), but it has closer ties to recent 
theoretical predictions from binary population synthesis models. 
For example, Ruiter et al. (2013) found that a “violent merger” 
DD model predicts at -1 power law shape for long delay SNe, 
but also includes a very prompt component that arises from a 
distinct subset of binary systems. A separate prompt channel 
for SN la explosions could also arise from a single degenerate 
pathway with a helium star donor (Wang et al. 2009; Claeys 
etal. 2014). 

8.3. DTD Model Fitting Results 

To find the most likely values for our two parameters ^i a and 
fp , we use three SN la rate data sets. First, we define the “All” 
data set, utilizing all available (non-redundant) volumetric rate 
measurements from the literature (see Graur et al. 2014, for 
a compilation table). Second, our “Ground” subsample picks 
out the 13 independent rate measurements from ground-based 
surveys. Finally, our “CANDELS+CLASH” sample isolates 
those two companion HST surveys. 

The first three columns of Table 6 summarize the maximum 
likelihood values for our DTD parameters rj i a and fp, when fitting 
to each of these subsamples. When using all of the available 
SN Ia survey data, we find ^ Ia = (1 .60 st ^o. ' 236 ^ 0 . 59 ) x SN Ia 

yr -1 M 0 ~ l and f P = 0.53 st ^;^ sy ^. This is a statistically 

acceptable fit, with a reduced x 2 °f 0.9 (p-value = 0.67). 

Fitting to the ground-based data alone, we find very similar 
best-fit parameters, with the prompt SN Ia fraction inching 
up to f P =0.59 and the efficiency remaining at rji a ~ 1.5. 
When we isolate the HST CANDELS and CLASH surveys, we 
get much larger uncertainties, but perhaps also a subtle hint 
at tension between the ground- and HST - based results: from 
the CANDELS+CLASH sample we get f P =0.21 st ^ s j^. 
The difference in these best-fit parameters reflects a (very) mild 
disagreement between the ground-based, primarily low-z rate 
measurements and the high-z constraints from HST. 

The source of this deviation is easily seen in Figure 11, 
where we plot two SNR(z) curves derived from the bifurcated 
DTD model. The (magenta) solid line shows the best fit to the 
ground-based data alone, with f P = 0.6. The (green) dashed 
line sets the prompt fraction to 20%, the best fit value for the 
CANDELS+CLASH data sample. These two HST surveys find 
a relatively low SN Ia rate at all redshifts z 1, which pulls 
the best-fit curve downward at high redshift, leading to the low 
best-fit fp. We can also see the slight tension between ground and 
HST measures in Figure 12, where we show confidence regions 


in the rji & vs. f P parameter space. The 68% contours from the 
ground- and ffi>T-based surveys fall just short of overlapping. 
This discrepancy is only slightly above la in significance, and 
comes with all the caveats cited above regarding the method 
for combining data from disparate surveys. Nevertheless, these 
HST data do sample the redshift range with the greatest leverage 
for constraining f P , so the scarcity of high-z SN Ia detections in 
multiple HST surveys should not be discounted. 

Table 6 and Figure 12 also present the total number of 
SNe Ia per stellar mass, Afi a / M * . This is computed by integrating 
the SN Ia rate over a Hubble time, and dividing by the total mass 
of formed stars. For the denominator, we take the integral of the 
Behroozi et al. (2013) CSFR(z) curve from Figure 10 (which 
assumes a Chabrier 2003 stellar initial mass function). To get the 
numerator — the total number of SN Ia explosions in a Hubble 
time — we can integrate the best-fit DTD-based SNR(z) model 
for each subsample of rate measurements. Those values are 
reported in the fourth column of Table 6. In the fifth column 
we list an alternative calculation, now directly integrating the 
SNR(z) data, without reference to any DTD model. The latter 
approach yields a much less precise constraint, but it is more 
appropriate for use as a test of progenitor models, because it does 
not presuppose any particular shape for the DTD. Note that we 
do not measure a data-only constraint from the ground-based 
subsample because it does not reach a high enough redshift. 

Figure 12 shows a color map in the background, reflecting 
the variation of Afi a /M* within the 17 i a — f P plane (assuming 
that the DTD follows our two-component toy model). The 
contours derived from both the ground-based and HST surveys 
are roughly aligned along lines of constant Afi a /M* (a single- 
color ridge in the color map). Hence the tight model-dependent 
constraints on Afi a /M * as reported in Column 4 of Table 6. 

All of the above 27i a /M* measurements are consistent with 
a value of roughly 1 x 10 _3 M 0 _1 . This is fully consistent 
with past measures of the volumetric rates, using similar stellar 
initial mass function (IMF) assumptions (e.g., Graur et al. 201 1). 
Other observational constraints, such as cluster SN Ia rates, have 
recently found values closer to 2 x 10 -3 M 0 _1 (Maoz & Badenes 
2010) — still consistent within the large error bars. However, 
theoretical predictions from binary population synthesis models 
are frequently lower by factors of 10 or more (Bours et al. 2013). 
This discrepancy between theory and observation remains one 
of the key concerns in the SN Ia progenitor problem. 

8.4. Interpretation and Speculation 

As described above, our analysis of all available SN Ia 
rate measurements suggests that the fraction of SN Ia explo- 
sions occurring <500 Myr after formation is fp ~ 50%. This 
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Figure 11. Comparing observed SN rates against DTD models. Gray markers show the collection of rates from the literature, using filled points for the ground-based 
surveys and open symbols for past HST surveys. As in Figure 9, large purple diamonds show the CLASH rates from Graur et al. (2014) and large orange squares show 
the CANDELS rates from this work. In this figure, vertical error bars show only the statistical uncertainties. Two curves show the SN la rates predicted by assuming a 
DTD that is proportional to t~ l for all times above 500 Myr, but with different assumptions for the fraction of SNe la that are prompt. The magenta solid line shows 
the best fit to the ground-based data (solid gray points), which has ~60% of all SNe la exploding within 500 Myr of birth. The green dashed line is the best-fit model 
when using the CANDELS +CL ASH data alone, for which the prompt SN la fraction is ~20%. 



_l I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


fp : prompt SN la fraction 

Figure 12. Constraints on the DTD normalization factor and the fraction of 
SNe la that are prompt explosions. Contours show the 68% and 95% confidence 
regions for the baseline assumptions (mid-dust, mid-rates) in the rj i a vs. fp 
parameter space. The background color map indicates the time-integrated 
SN Ia efficiency, N\ a /M *, for each point in that parameter space. Dashed 
contours show the confidence regions derived from only HSTdata, reaching 
to z ~ 2.5. Solid contours are from the collection of all ground-based SN 
surveys, dominated by measurements at z < L 

observed value of fp is consistent with the simplistic assump- 
tion of a t -1 DTD that continues without truncation all the way 
down to 40 Myr, which yields f P = 0.43. A prompt fraction 
close to 50% is also observationally supported by several lines 


of evidence in the local universe. Mannucci et al. (2006) first 
proposed that roughly half of all SN Ia explode promptly after 
formation, based on observations of SN Ia host galaxies at low 
redshift. Building on that work, Raskin et al. (2009) used mea- 
surements of low-z SN Ia environments on sub-galactic scales to 
infer that most of those prompt SNe Ia explode at 200-500 Myr. 
Mennekens et al. (2013) used binary population synthesis (BPS) 
to predict the distribution of chemical enrichment in our galaxy 
over time. Comparing this to observations of [Fe/H] in nearby 
G dwarfs, they infer that prompt explosions must make up a 
large fraction of the SN Ia population (and thereby contribute to 
rapid galactic enrichment). 

Can this measurement of the prompt SN Ia fraction be used 
to distinguish SD and DD progenitor models? BPS calculations 
generally agree that SD pathways preferentially generate prompt 
SN Ia explosions. In particular, SD progenitor models in 
which the companion is a naked He star are found to peak 
at t ~ 100 Myr after formation, while those with a normal 
main sequence or giant companion preferentially explode at 
200-500 Myr (Wang et al. 2009; Mennekens et al. 2010; 
Greggio 2010; Claeys et al. 2014). Some BPS modeling also 
finds that DD progenitors could contribute substantially to the 
population of SN Ia explosions younger than 500 Myr (Ruiter 
et al. 2009; Greggio 2010; Ruiter et al. 2013), although recent 
work by Claeys et al. (2014) suggests that the DD pathway 
does not dominate the DTD until t > 500 Myr. At the 
moment, we can only say that a prompt fraction fp ~ 50% 
is commonly predicted by models that include both SD and DD 
progenitors — but it may be possible in a pure DD model as well. 

When we analyze the HST sample in isolation, we find a hint 
that f P could be closer to ~20%, though further analysis and a 
larger SN Ia sample will be needed to improve this measurement. 
Let us indulge in a bit of speculation and suppose that these 
future improvements confirm that/p is close to 20% at high 
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redshifts. We would then need a theoretical explanation for 
how the prompt component of the SN la population could be 
suppressed at redshifts z > 1.5. One possible explanation would 
be that the prompt SN la component is dominated by an SD 
progenitor pathway with a strong metallicity dependence. Such 
a metallicity dependence has been proposed by requiring an 
optically thick wind from the WD to regulate the mass transfer 
rate, allowing the WD to grow to the Chandrasekhar mass 
limit (Hachisu et al. 1996, 1999). This wind would be absent 
in low metallicity WDs, resulting in a suppression of prompt 
SN la explosions in the early universe (Kobayashi et al. 1998; 
Kobayashi & Nomoto 2009). No such metallicity threshold 
is expected for DD models, so this scenario could provide a 
clean way to disentangle the relative contributions of those two 
progenitor pathways. 

9. SUMMARY 

We have presented a sample of 65 SNe from the five 
CANDELS fields. This sample, collected in concert with the 
CLASH SN search, is the first to extend SN la detections beyond 
z = 2, and the first to detect SNe at z ~ 1.5 in rest- frame 
optical bands. These SNe have been classified primarily through 
the application of STARDUST, a new Bayesian photometric 
classifier that is optimized for working with light curves of high 
redshift SNe. We have spectroscopic redshifts from the SN and/ 
or host galaxy for 82% of the sample (53 of 65), although we 
rely on photometric redshifts for 43% of the SNe at z > 1.5 
(6 of 14). Our SN classification probabilities are in general 
tightly constrained by well-sampled light curves, rest-frame 
UV-optical colors, and well-defined redshifts. The primary 
sources for potential systematic biases in our classifications are 
(1) a redshift-dependent prior describing the relative fraction 
of SNe that are Type la, and (2) the assumed distribution 
of dust extinction values. For the former, a test using host 
galaxy information to replace the class prior indicates that our 
systematic uncertainty estimates are appropriate. 

From the CANDELS SN sample we have measured the 
volumetric SN la rate in five redshift bins reaching to z = 
2.5. We find that the CANDELS SN la rate measurement at 
z ~ 1.25 is a factor of 2 lower (~2 sigma) than past HST rate 
measurements at the same redshift (Dahlen et al. 2008; Barbary 
et al. 2012), but is consistent with the concurrently measured 
SN la rates from CLASH (Graur et al. 2014). We attribute this 
discrepancy to Poisson noise, due to the very small sample sizes 
in all of these HST surveys. At higher redshifts the CANDELS 
rate measurements remain approximately flat with redshift. 

Combining these CANDELS results with other surveys from 
HST and from the ground, we have examined the constraints 
that can be placed on SN la DTD models (Table 6, Figure 11). 
We have invoked a simple two-component model with a t~ l 
distribution for long delay times ( t > 500 Myr), and a constant 
rate at shorter times. We find that the ground-based rates 
(primarily at z < 1) and the full sample of all available 
SN la rates are both best matched with a prompt SN la fraction of 
fp = 0.5. When the CANDELS+CLASH surveys are analyzed 
in isolation, the best fit for the prompt SN la fraction falls to 
fp = 0.21^^. This is much lower, but with very large 
error bars. 

Collectively, the constraints from all available volumetric 
SN data indicate that the prompt SN la fraction cannot be 
much larger than 60%. As described above, there is room for 
substantial improvement in our measurement of the prompt 
SN la fraction — without needing to acquire more data. System- 


atic uncertainties can be reduced (or at least better understood) 
by combining the existing HST surveys in a composite analysis 
that handles SN classification and dust obscuration in a con- 
sistent manner. Increasing the high-z SN la sample size would 
also help, of course. The HST Frontier Fields initiative will pro- 
vide the next opportunity for new high-z SN discoveries. This 
program will utilize ~900 orbits of HST observations for very 
deep imaging of six massive galaxy clusters over three years. 
An approved HST program for SN discovery and follow-up 
(PI: Rodney, PID: 13386) is expected to deliver an SN la sam- 
ple that reaches out to z ~ 3. 

This work is dedicated to the memory of our dear friend and 
colleague, Tomas Dahlen. 

The authors would like to thank the referee, Massimo Della 
Valle, for helpful comments and suggestions. We also thank 
the STScI science support staff for their extraordinary efforts 
in executing the MCT programs, especially program coordi- 
nators Patricia Royle and Beth Perriello, as well as the entire 
STScI scheduling team that made our HST ToO program possi- 
ble. Thanks also to Ori Fox, Patrick Kelly, Isaac Shivvers, Brad 
Tucker, and WeiKang Zheng for assistance with Keck observa- 
tions; to Ismael Botti, Alice Mortlock, and Omar Almaini for 
redshift and AGN classifications of SN candidates in the UDS 
field; to Stephane Blondin for VLT spectral reductions and as- 
sistance with SNID modifications; to Jennifer Lotz and Ryan 
Foley for discussions improving the measurement and applica- 
tion of SN host galaxy morphologies; and to Thomas Holoien 
and Mark Ziegler for “fake SN” searching assistance. 

The CANDELS program was supported by NASA through 
grant HST-GO- 12060, and the SN-MCT program through grant 
HST-GO- 12099 from the Space Telescope Science Institute, 
which is operated by Associated Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. Support 
for S.A.R. was provided by NASA through Hubble Fellowship 
grant HST-HF-51312. Support for this research at Rutgers 
University was provided in part by NSF CAREER award AST- 
0847157 to S.W.J. The Dark Cosmology Centre is supported 
by the Danish National Research Foundation. R.P.K. thanks 
the National Science Foundation for AST-1211196, and the 
John Simon Guggenheim Foundation for support. This work 
was supported by NASA Keck PI Data Awards (to Rutgers 
University, PI: S.W.J.), administered by the NASA Exoplanet 
Science Institute. J.M.S. is supported by an NSF Astronomy 
and Astrophysics Postdoctoral Fellowship under award AST- 
1302771. A.V.F. is also grateful for the support of the National 
Science Foundation (NSF) grant AST-121 1916, the TABASGO 
Foundation, and the Christopher R. Redlich Fund. 

This research was based primarily on observations made 
with the NASA/ESA Hubble Space Telescope , delivered by 
the data archive team at the Space Telescope Science Institute 
(STScI), which is operated by the association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 
5-26555. We used the HST Science Archive hosted by the 
Canadian Astronomy Data Centre (CADC/NRC/CSA), as well 
as the Mikulski Archive for Space Telescopes (MAST). 

This work is based in part on observations made with ESO 
telescopes at the La Silla Paranal Observatory under program 
IDs 086.A-0660 and 088.A-0708. 

Some of the data presented herein were obtained at the 
W. M. Keck Observatory from telescope time allocated to NASA 
through the agency’s scientific partnership with the California 
Institute of Technology and the University of California. The 


18 


The Astronomical Journal, 148:13 (28pp), 2014 July 


Rodney et al. 


Table 7 

Host Galaxy Likelihood Distributions for Defining 
Host-based Classification Priors a 


Category 

P(D|Ia) 

P(D |CC) 

Morphology 

spheroid (E,S0) 

0.25 ± 0.03 

0.01 ± 0.01 

spheroid+disk (S0,Sa) 

0.26 ± 0.03 

0.15 ±0.02 

disk (Sb,Sbc,Sc) 

0.35 ± 0.03 

0.53 ± 0.03 

disk+irregular (Sc,Scd) 

0.10 ±0.02 

0.20 ± 0.02 

irregular (Scd,Irr) 

0.04 ± 0.01 

0.11 ±0.02 

SED Type b 

Passive (B - K > 3.75) 

0.35 ± 0.06 

0.08 ± 0.06 

Active (2.75 < B - K < 3.75) 

0.52 ± 0.08 

0.62 ± 0.08 

Starburst (B — K < 2.75) 

0.13 ±0.04 

0.30 ± 0.04 


Notes. 

a Likelihood estimates follow the galsnid for malism (Foley & Mandel 2013), 
using LOSS host galaxy data compiled in Leaman et al. (201 1). 
b Approximate rest-frame B — K colors are given in Vega mags. In AB mags the 
thresholds are 0.82 and 1.82. 

Observatory was made possible by the generous financial 
support of the W. M. Keck Foundation. The authors wish to 
recognize and acknowledge the very significant cultural role 
and reverence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. We are most 
fortunate to have the opportunity to conduct observations from 
this mountain. 

This work is based in part on observations obtained at the 
Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini part- 
nership: the National Science Foundation (United States), the 
National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), Ministerio da Ciencia, 
Tecnologia e Inova^ao (Brazil) and Ministerio de Ciencia, Tec- 
nologia e Innovation Productiva (Argentina). Contributing data 
came from Gemini programs GN-2011A-Q-14, GN-2011B-Q- 
18, GN-2012A-Q-32, GN-2013A-Q-25, GS-2011A-Q-16, GS- 
2011B-Q-18, GS-2012A-Q-17, and GS-2013A-Q-19. 

This research has made use of the VO Datascope, developed 
with the support of the National Science Foundation under 
Cooperative Agreement ASTO 122449 with the Johns Hopkins 
University, and hosted by the Astrophysics Science Division and 
the High Energy Astrophysics Science Archive Research Center 
(HE AS ARC), a service of Goddard Space Flight Center and the 
Smithsonian Astrophysical Observatory. This research has made 
use of the SIMBAD database, operated at CDS, Strasbourg, 
France. This research has made use of NASA’s Astrophysics 
Data System Bibliographic Services. This research made use 
of Astropy, a community-developed core Python package for 
Astronomy (Astropy Collaboration et al. 2013). 

Facilities: HST (WFC3), Keck (DEIMOS), Keck (LRIS), 
VLT (X-shooter), VLT (FORS2), Gemini (GMOS-N), Gemini 
(GMOS-S) 

APPENDIX A 

THE HOST-BASED CLASS PRIOR 

In this Appendix we perform a test to see if our systematic 
uncertainty estimates are accurately reflecting the bias that could 
arise from using an incorrect redshift-dependent prior P(Ia, z) in 
the STARDUST classifier. For this test we measure the change in 


Table 8 

Change in Observed SN la Counts and the Volumetric Rate 
When Adopting the Host Galaxy Prior 


Redshift 

A/Vobs 

ASNR 

% of Sys. Err. 

0.25 

+0.13 

+0.03 

26 

0.75 

+1.99 

+0.14 

62 

1.25 

+0.93 

+0.07 

21 

1.75 

+0.41 

+0.05 

10 

2.25 

+0.33 

+0.13 

28 


the volumetric SN la rate that occurs when we adopt a redshift- 
independent prior based on host galaxy data. 

It is well established that CCSNe are almost never observed 
in galaxies dominated by older stellar populations (passive, 
red ellipticals). The rate of SNe la per unit mass increases in 
galaxies with young stellar populations, but not as sharply as the 
specific CCSN rate. The galsnid SN classifier (Foley & Mandel 
2013) exploits these relationships between SNe and their host 
environments to define a posterior classification probability that 
is completely independent of the SN photometry. 

The galsnid approach relies on a database of observed SN 
host galaxy properties (morphology, color, luminosity, etc.) 
from the Lick Observatory Supernova Survey (LOSS; Leaman 
et al. 2011). These data are used to define a set of likelihoods 
P(T>h | Type) giving the probability of observing a set of host 
galaxy properties, D^, if the SN hosted in that galaxy is of the 
given type. The galsnid posterior classification probability for 
Type la is 

n 

P(Ia|D ft ) = ^- 1 J P(Ia) P(Di |Ia), (Al) 

1 = 1 

where P(Ia) is the la classification prior, Dhj are the n observed 
host galaxy properties, and kh is a normalization term, which 
ensures that the sum of posterior classification probabilities for 
all SN la and CCSN types is unity. For our implementation 
of galsnid , the prior P (la) — and the corresponding priors for 
CCSNe — are fixed to match the observed fraction of SNe in 
a volume-limited sample that are of each type: P (la) = 0.25, 
P (Ib/c) = 0.19, and P (II) = 0.57 (Smartt et al. 2009; Li 
et al. 2011). For the host galaxy observables, D^, we employ 
the two quantities that provide the strongest discriminatory 
power according to Foley & Mandel (2013): morphology and 
color. 

Table 7 translates from the CANDELS categories for host 
galaxy morphology (spheroid, disk, or irregular) into their 
approximate counterparts on the Hubble sequence (E, SO, etc.). 
That table also describes the translation from the CANDELS 
SED types (passive, active, or starburst) into the rest-frame 
( B — K) color. 

For SNe in faint hosts where the host galaxy morphology is 
undefined, the galsnid posterior reflects only the constraint from 
color/SED type. In the two cases where the SN host galaxy is 
totally undetected, the photometric redshifts are unconstrained 
and galsnid simply reflects back the input prior, P (la) in 
Equation (Al). The likelihood values in Columns 2 and 3 of 
Table 7 mimic Table 1 of Foley & Mandel (2013) and provide 
all the information needed to determine P (la, host) for any 
CANDELS SN host galaxy. 

A.l. Combining galsnid and STARDUST 

To incorporate the galsnid information into STARDUST, we 
adopt the galsnid posterior, F(Ia|Dh), as a redshift-independent 
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Table 9 

51 Supernovae with z < 1.5 


Name 

R.A. 

Decl. 

P (Ia|ZU) a 

P (Ia|£>h 0st ) b 

^SN C 

(±) 

z Source d 


(J2000) 

(J2000) 






GSDIOTum 

03:32:17.705 

-27:50:57.50 

0 00 +0 0 ° 
uuu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

0.124 

(0.001) 

host spec-z 

COS12Cli 

10:00:15.958 

+02:12:36.27 

n 29 +014 
-0.28 

0 29 +0 - 03 

-0.29 

0.187 

(0.001) 

host spec-z 

GSD12Roc 

03:32:06.368 

-27:47:26.63 

0 00 +0 0 ° 
u - uu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

0.346 

(0.003) 

host spec-z 

EGSllTho 

14:19:31.775 

+52:51:56.16 

f) QQ +0.01 
uyy -0.98 

0 QQ + 0 00 
-0.98 

0.354 

(0.001) 

host spec-z 

COS12Ken 

10:00:36.010 

+02:15:26.53 

0.00^“ 

0 00 +0 0 ° 
u - uu -0.00 

0.373 

(0.001) 

host spec-z 

GND13She 

12:36:09.877 

+62:14:06.00 

0+8 "oS 

0 7 1 + 0 - 25 
U - :)A -0.31 

0.473 

(0.001) 

host spec-z 

GSDllTay 

03:32:54.505 

-27:47:04.23 

0 67 +0 - 21 
U - 0/ -0.64 

n 70 +0.12 
U - /Z -0.64 

0.535 

(0.001) 

host spec-z 

GSDllFor 

03:32:14.300 

-27:47:13.30 

n 41 + 0 - 26 
U ^ A -0.41 

0 47 + 013 

U-+-5 _ 0 42 

0.578 

(0.001) 

host spec-z 

GND13Bid 

12:36:41.325 

+62:11:42.21 

0 00 +0 0 ° 
uuu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

0.585 

(0.001) 

host spec-z 

COS12Eis 

10:00:47.233 

+02:11:50.50 

f) ci +0.14 
_ 0 . 18 

0-69 IT 

0.605 

(0.001) 

host spec-z 

UDSllGar 

02:17:39.633 

-05:11:37.14 

0 03 +0 09 

u - u:> - 0.03 

0-08 Tos 

0.651 

(0.001) 

host spec-z 

GSDllWor 

03:32:10.739 

-27:48:07.17 

1 00 +00 ° 
AUU -0.00 

1 00 +00 ° 
A - uu -0.00 

0.654 

(0.001) 

host spec-z 

GSDllRoo 

03:32:31.587 

-27:46:12.53 

0 00 +0 0 ° 
uuu -0.00 

n ni +0.00 

u.ui _ 000 

0.655 

(0.001) 

host spec-z 

GSWllJac 

03:32:32.300 

-27:54:20.55 

0 17 +013 
U - A ' -0.11 

n 1 0 + 0.04 
- 0.06 

0.659 

(0.001) 

host spec-z 

COS12Tru 

10:00:38.322 

+02:11:36.41 

o 46 + 0-32 
u -^° - 0.46 

0 52 +0 - 22 

U-JZ _ 0 52 

0.665 

(0.001) 

host spec-z 

COS12Rea 

10:00:31.917 

+02:14:16.15 

0 00 +0 0 ° 
u - uu —0.00 

0 00 +0 0 ° 
u.uu _ 0 00 

0.679 

(0.001) 

host spec-z 

GSDllBuc 

03:32:28.800 

-27:52:32.00 

0.01 To 1 ! 

0 1 n +006 
U - AU -0.10 

0.679 

(0.001) 

host spec-z 

GSDllHar 

03:32:30.441 

-27:45:18.70 

0-32 T 3 2 

0 25 +012 

U - ZJ - 0.25 

0.681 

(0.001) 

host spec-z 

COS 12 Aid 

10:00:15.252 

+02:17:32.12 

0 00 +0 0 ° 
U-UU -0.00 

0 00 +0 0 ° 
u.uu _ 0 00 

0.731 

(0.001) 

host spec-z 

GSDllLin 

03:32:29.812 

-27:49:19.71 

0 00 +0 0 ° 
u - uu -0.00 

n oi +0 00 
U - UA -0.01 

0.734 

(0.001) 

host spec-z 

GSDllAda 

03:32:19.805 

-27:54:10.04 

n f)1 + 002 
UUA -0.01 

0 01 +001 
U - UA -0.01 

0.735 

(0.001) 

host spec-z 

COS12Mon 

10:00:26.737 

+02:15:13.74 

0 07 +0 06 
u - uz _ 0 02 

0 4Q + 019 
U-+^ —O.49 

0.777 

(0.001) 

host spec-z 

GND13Ful 

12:36:19.201 

+62:15:12.58 

0 00 +0 0 ° 
uuu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

0.783 

(0.001) 

host spec-z 

GND12Bur 

12:36:32.536 

+62:15:32.62 

0 00 +0 0 ° 
u - uu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

0.82 

(0.13) 

host+SN phot-z 

GND12Daw 

12:36:41.340 

+62:18:52.50 

0 00 +0 0 ° 
uuu _o.oo 

0 00 +0 0 ° 
u.uu _ 0 00 

0.83 

(0.02) 

SN spec-z + phot-z 

GND13Vin 

12:37:06.354 

+62:15:17.79 

0 1 1 +018 
UAA -0.11 

0 1 3 + 009 

U - A:> - 0.13 

0.840 

(0.001) 

host spec-z 

GNW12Ger 

12:37:17.023 

+62:20:38.67 

n nc +0.02 
U-U- 3 - 0.03 

0 66 +0 0 ° 
u -°° -0.10 

0.850 

(0.001) 

host spec-z 

GND13Roy 

12:36:46.223 

+62:15:27.13 

0 00 +0 0 ° 
uuu -0.00 

0 00 +0 0 ° 
u - uu -0.00 

0.851 

(0.001) 

host spec-z 

GSDllVan 

03:32:19.037 

-27:47:17.90 

0 05 +012 

U ‘ UJ - 0.05 

0 19 +018 

u - 0.19 

0.886 

(0.097) 

host+SN phot-z 

EGSllPer 

14:18:28.318 

+52:42:45.88 

0 QQ + 0 00 
u.vv _ 0 01 

1 00 +00 ° 
A - uu -0.00 

0.915 

(0.001) 

host spec-z 

COS12Mik 

10:00:31.687 

+02:26:13.73 

0 00 +0 0 ° 
uuu -0.00 

0 01 +00 ° 
U - UA -0.00 

0.927 

(0.001) 

host spec-z 

GND12Cal 

12:37:10.487 

+62:15:47.67 

1 00 +00 ° 
AUU -0.00 

1 00 +00 ° 
A - uu -0.00 

0.941 

(0.005) 

host spec-z 

GND13Wol 

12:36:58.946 

+62:18:10.15 

n 2 1 + 0.08 
UJA -0.12 

0 +0.02 

u -7o _o.05 

0.943 

(0.001) 

host spec-z 

GSDllMad 

03:32:18.781 

-27:52:42.05 

1 00 +00 ° 
AUU -0.00 

1 00 +0 - 00 
A - uu -0.00 

0.988 

(0.001) 

host spec-z 

GNW12Gor 

12:36:20.704 

+62:08:45.08 

n oq +01 ° 
uuy - 0.08 

0 23 +0 - 07 

u - z:> - 0.15 

1.016 

(0.001) 

host spec-z 

GND13Reh 

12:36:54.452 

+62:11:52.47 

0 00 +0 0 ° 
uuu -0.00 

0 09 +0 03 
u>uy - 0.06 

1.019 

(0.001) 

host spec-z 

GND13Jay 

12:36:41.380 

+62:11:30.02 

1 00 +00 ° 
AUU -0.01 

1 00 +00 ° 
A - uu -0.00 

1.03 

(0.01) 

host spec-z 

GND13Gar 

12:36:40.806 

+62:11:14.16 

1 00 +00 ° 
1 uu -0.00 

1 00 +00 ° 
A - uu -0.00 

1.067 

(0.008) 

SN spec-z + phot-z 

COS12Her 

10:00:47.446 

+02:15:17.74 

o.oi T& 

0.02 T01 

1.082 

(0.077) 

host+SN phot-z 

GSD12Agn 

03:32:25.902 

-27:50:19.62 

0.80 T 3 ! 

0.82 T 7 3 

1.095 

(0.001) 

host spec-z 

EGS13Tan 

14:20:02.098 

+53:00:16.36 

0+5 Tos 

0.35 Tt 

1.137 

(0.439) 

host+SN phot-z 

GND12Tom 

12:36:42.543 

+62:18:21.93 

0.86 T 9 5 

0.85 

1.14 

(0.001) 

host spec-z 

EGS13Wai 

14:20:28.534 

+53:04:58.61 

1 00 +000 
AUU -0.00 

1 00 +00 ° 
A - uu -0.00 

1.166 

(0.001) 

host spec-z 

GND13Cam 

12:37:07.354 

+62:10:26.90 

1 00 +00 ° 
A-UU _ 001 

1 00 +00 ° 
A-UU _ 000 

1.222 

(0.002) 

host spec-z 

GSDllWas 

03:32:20.856 

-27:49:41.48 

1 00 +00 ° 
AUU -0.00 

1 00 +00 ° 
A - uu -0.00 

1.33 

(0.02) 

SN spec-z + phot-z 

GND12Whe 

12:36:24.220 

+62:17:39.70 

0 00 +0 0 ° 
u - uu -0.00 

0-00 1° 0 % 

1.335 

(0.005) 

host spec-z 

EGS13Cha 

14:20:24.200 

+53:03:37.50 

n 26 +0 - 01 

U ~ 3 ° _q 14 

0.80 TA 

1.40 

(0.01) 

host spec-z 

EGSllOba 

14:20:32.663 

+53:02:48.18 

0 Q 1 +0 02 
U ^ A - 0.17 

0.90 

1.409 

(0.001) 

host spec-z 

EGSllNix 

14:20:48.603 

+53:00:26.47 

n 40 + 0 01 
u.+u _ 0 n 

0.58 !°„% 

1.451 

(0.077) 

host+SN phot-z 

GND12Kin 

12:37:13.005 

+62:16:30.83 

0 00 +0 0 ° 
u.uu _ 0 00 

0 00 +0 0 ° 
U.UU _ 0 00 

1.46 

(0.07) 

host+SN phot-z 

GND13Fin 

12:37:18.173 

+62:17:09.65 

0 07 +0 06 
u.uj _ 0 03 

0 07 +0 01 

U.UZ _ 0 01 

1.488 

(0.001) 

host spec-z 


Notes. 

a Type la SN classification probability from STARDUST, using the redshift-dependent class prior. Uncertainties reflect systematic biases due to the class prior and extinction 
assumptions (Sections 4.2 and 4.3). 

b Type la SN classification probability from STARDUST, using the galsnid host galaxy prior. Uncertainties reflect systematic biases due to the class prior and extinction assumptions. 
c Posterior redshift and uncertainty, as determined by the STARDUST light curve fit. 

d The host/SN values indicate whether the redshift is derived from the host galaxy, the SN itself, or a combination; spec-z/phot-z specify a spectroscopic or photometric redshift. A 
value of host+SN phot-z means the redshift is derived from a STARDUST light curve fit, with the host galaxy phot-z used as a prior. 
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GSDIOTum, z = 0.1239 COS12CH, z = 0.187 GSD12Roc, z = 0.43 ± 0.22 



EGSllTho, z = 0.354 COS12Ken, z = 0.371 GND13She, z = 0.473 





GSDllWor, z= 0.654 


COS12Eis,z= 0.605 


UDSllGar,z= 0.651 



COS12Rea, z= 0.679 GSDllBuc, z= 0.6787 GSDllHar, z= 0.681 




GND12Daw, z = 0.83 ± 0.02 GND13Vin, z = 0.8398 GNW12Ger, z = 0.850 



Figure 13. Detection images for 27 SN from the CANDELS fields with redshifts z 0.85. Each image triplet shows H band (F160W) images with the template 
image on the left, the discovery epoch image in the middle and the difference image on the right. All images have a width of about 6 arcsec, with north up and east to 
the left. The position of the SN is marked by (red) crosshairs in every frame. 

(A color version of this figure is available in the online journal.) 
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EGS13Wai, z= 1.1662 GND13Cam, z= 1.222 GSD12Was, z = 1.33 ± 0.02 




Figure 14. Detection images for 24 SN from the CANDELS fields with redshifts 0.85 < z < 1.5. Each image triplet shows //band (F160W) images with the template 
image on the left, the discovery epoch image in the middle and the difference image on the right. All images have a width of about 6 arcsec, with north up and east to 
the left. The position of the SN is marked by (red) crosshairs in every frame. 

(A color version of this figure is available in the online journal.) 


prior P (la, host), replacing P( la, z) in Equation (2). Note 
that this assumption is assuredly incorrect: at redshift z ~ 2 
the fraction of SNe la appearing in red early-type galaxies 
must be much lower than it is locally, simply because there 
are far fewer of those passive old galaxies at high redshift. 
However, for the purpose of this systematic test, we will 
make the brazen assumption that the relationships between 


SN types and their host galaxy properties do not evolve with 
redshift. 

To define our baseline rate measurement, we have relied on 
the redshift-dependent mid-rate P( la, z) prior (the green solid 
curve in Figure 6). Because the galsnid - based prior P (la, host) 
prior is redshift independent, we can use it to check for strong 
redshift biases in the rate-based P( la, z) prior. In principle, these 
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GSDIOTum 

* F350LP^F160WD F125W 



P(Ia|D)=0.00 
P(Ibc|D)=0.14 
P(II|D)=0.86 
z sn =0.124± 0.000 

Best-fit Model: 

IIP(SDSS-018297) 

MJD pk =55535 

z=0.124 

A v =1.00 

R v =3.1 

Am=0.7 

xH-7 


COS12CH 

Of160wD F125W> F814W^ F606W 



P(Ia | D)=0.29 
P(Ibc|D)=0.00 
P(II|D)=0.71 
z SN =0.187± 0.001 

Best-fit Model: 

IIP(SDSS-015031) 

MJD pk =55948 

z=0.187 

A v =2.66 

R v =3.1 

Am=0.9 

xH-8 


GSD12Roc 

♦ F160W A F850LP* F350LP 



P(Ia | D)=0.00 
P(Ibc | D)=0.00 
P(II | D)=l .00 
z sn = 0-346± 0.003 

Best-fit Model: 

IIP(SDSS-015031) 

MJD pk =55979 

z=0.346 

A v =0.48 

R v =3.1 

Am=0.6 

xH-° 


EGSllTho 

Of160wB F125W^ F814W^ F606W 


P(Ia|D)=0.99 

P(Ibc|D)=0.00 

P(II|D)=0.01 



COS12Ken 

♦ F160WD F125Wl> F814W^ F606W 



P(Ia | D)=0.00 
P(Ibc|D)=0.01 
P(II|D)=0.99 
z SN =0.371± 0.001 

Best-fit Model: 

IIP(SDSS-015320) 

MJD pk =55961 

z=0.371 

A v =0.00 

R v =3.1 

Am=0.6 

Xv=30 


GND13She 

Of160WD F125W* F350LP 



P(Ia|D)=0.18 
P(Ibc|D)=0.81 
P(II|D)=0.00 
z sn =0.473± 0.001 

Best-fit Model: 

Ic(SDSS-004012) 

MJD pk =56327 

z=0.473 

A v =0.00 

R v =3.1 

Am=0.9 

Xv=2-0 


GSDllTay 

Of160W« F125W- F606W* F350LP 



P(Ia|D)=0.67 
P(Ibc|D)=0.16 
P(II | D)=0. 1 8 
z SN =0.535± 0.001 

Best-fit Model: 

Ia(SALT2) 

MJD pk =55839 

z=0.535 

xi=-0.16 

C=0.5 

(3=4.1 

Xv=0-4 


GSDllFor 

Of160WA F850LP- F350LP 



P(Ia|D)=0.41 
P(Ibc|D)=0.37 
P(II|D)=0.22 
z sn =0.578± 0.001 

Best-fit Model: 

Ic(SDSS-015475) 

MJD pk =55876 

z=0.578 

A v =0.97 

R v =3.1 

Am=0.3 

xH-7 


GND13Bid 

♦ f160WD F125W- F814W- F350LP 



P(Ia | D)=0.00 
P(Ibc|D)=1.00 
P(II|D)=0.00 
z sn =0.585± 0.001 

Best-fit Model: 

Ic(CSP-2004gq) 

MJD pk =56329 

z=0.585 

A v =2.90 

R v =3.1 

A m =-0.7 

xH-6 


COS12Eis 

OF160WH F125W^ F814W- F606W 



P(Ia|D)=0.53 
P(Ibc|D)=0.14 
P(II|D)=0.33 
z sn =0.605± 0.001 

Best-fit Model: 

Ia(SALT2) 

MJD pk =55980 

z=0.605 

X!=-0.79 

C=-0.1 

(3=4.1 

Xv=0.6 


UDSlOGar 

Of160wH F125wl> F814W^ F606W 


P(Ia | D)=0.03 
P(Ibc|D)=0.09 
P(II|D)=0.87 



GSDllWor 

Of160WA F850LP^ F606W 



P(Ia|D)=1.00 
P(Ibc | D)=0.00 
P(II|D)=0.00 
z sn =0.654± 0.001 

Best-fit Model: 

Ia(SALT2) 

MJD pk =55550 

z=0.654 

X!=-0.47 

C=0.1 

(3=4.1 

xH-6 


GSDllRoo 

Of160W- F850LP® F350LP 



P(Ia|D)=0.00 
P(Ibc|D)=0.02 
P(II|D)=0.98 
z sn =0.655± 0.001 

Best-fit Model: 

IIN(SDSS-012842) 

MJD pk =55799 

z=0.655 

A v =0.00 

R v =3.1 

Am=0.2 

xH-8 


GSWllJac 

Of160wD F125W* F350LP 



P(Ia|D)=0.17 
P(Ibc|D)=0.52 
P(II|D)=0.32 
z sn =0.659± 0.001 

Best-fit Model: 

Ic(CSP-2004fe) 

MJD pk =55605 

z=0.659 

A v =0.00 

R v =3.1 

Am=0.3 

X 2 V=0.0 


COS12Tru 

Of160WD F125W^ F814W < 



P(Ia | D)=0.46 
P(Ibc | D)=0.28 
P(II|D)=0.26 
z sn =0.665± 0.001 

Best-fit Model: 

Ib(CSP-2007Y) 

MJD pk =55947 

z=0.665 

A v =0.48 

R v =3.1 

Am=0.4 

X5=E5 


COS12Rea 

Of160WD F125W^ F814W^ F606W 



P(Ia|D)=0.00 
P(Ibc|D)=0.14 
P(II|D)=0.86 
z SN =0.679± 0.001 

Best-fit Model: 

IIP(SDSS-003818) 

MJD pk =55934 

z=0.679 

A v =0.48 

R v =3.1 

Am=0. 1 

XH.2 


GSDllBuc 

^F160WO F125W A F850LP^ F814W 



P(Ia|D)=0.01 
P(Ibc|D)=0.16 
P(II|D)=0.83 
z sn =0.679± 0.001 

Best-fit Model: 

IIP(SDSS-018441) 

MJD pk =55750 

z=0.679 

A v =0.24 

R v =3.1 

Am=-0.6 

xH-5 


GSDllHar 

<>F160Wn F125W- F350LP 


P(Ia|D)=0.32 

P(Ibc|D)=0.26 

P(II|D)=0.43 



Figure 15. STARDUST light curve matches for the first 18 SN from the CANDELS fields in redshift order, with redshifts z < 0.7, as in Figure 5. 

(A color version of this figure is available in the online journal.) 
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COS 12 Aid 

Of160WD F125W^ F814W* F606W 


P(Ia|D)=0.00 
P(Ibc | D)=0.05 
P(II|D)=0.95 



GSDllLin 

Of160WD F125W* F350LP 



P(Ia|D)=0.01 
P(Ibc|D)=0.24 
P(II|D)=0.75 
z sn =0.734± 0.001 

Best-fit Model: 

IIP(SDS S-0000 1 8) 

MJD pk =55814 

z=0.734 

A v =0.00 

R v =3.1 

Am=0. 1 

xH-0 


GSDllAda 

#F160W“ F850LP “ F350LP 



P(Ia|D)=0.01 
P(Ibc|D)=0.75 
P(II|D)=0.24 
z sn =0.735± 0.001 

Best-fit Model: 

Ib(SDSS-0 14492) 

MJD pk =55831 

z=0.735 

A v =0.48 

R v =3.1 

Am=-0.1 

xH-6 


COS12Mon 

Of160WD F125W^ F814W- F606W 



P(Ia|D)=0.02 
P(Ibc|D)=0.14 
P(II|D)=0.84 
z SN =0.777± 0.001 

Best-fit Model: 

IIP(SDSS-0 18590) 

MJD pk =55957 

z=0.777 

A v =0.00 

R v =3.1 

Am=0.1 

xH-i 


GND13Ful 

Of160WA F850LP® F350LP 



P(Ia|D)=0.00 
P(Ibc|D)=0.05 
P(II|D)=0.95 
z SN =0.783± 0.001 

Best-fit Model: 

IIN(SDSS-012842) 

MJD pk =56430 

z=0.782 

A v =0.48 

R v =3.1 

Am=-0.2 

Xv = F4 


GND12Bur 

O F 1 60W “ F850LP • F350LP 
□ F125W - F814W 




Best-fit Model: 

IIN(SDSS-013449) 

MJD„ k =56244 


56220 56250 56280 56310 56340 


GND12Daw 

Of160WO F125W^ F814W® F350LP 



P(Ia|D)=0.00 
P(Ibc | D)=0.00 
P(II| D)=l .00 

z sn =0.810± 0.008 

Best-fit Model: 

IIP(SDSS-000018) 

MJD pk =56103 

z=0.810 

A v =0.50 

R v =3.1 

Am =-1 .9 

xH-2 


GND13Vin 

Of160WE1 F125W® F350LP 
OF140W-F814W 




Best-fit Model: 

IIP(SDSS-000018) 

MJD ok =56448 


GNW12Ger 

O F140WO F105W- F814W* F350LP 


56430 56460 56490 56520 56550 



P(Ia|D)=0.05 
P(Ibc|D)=0.95 
P(II|D)=0.00 
z sn =0.850± 0.001 

Best-fit Model: 

Ic(SNLS-04D4jv) 

MJD pk =56175 

z=0.850 

A v =0.00 

R v =3.1 

Am=-0. 1 

xH-2 


GND13Roy 

<>F160WlI F125W^ F814W* F350LP 


P(Ia|D)=0.00 

P(Ibc|D)=0.03 

P(II|D)=0.97 



GSDllVan 

O F160W A F850LP ◄ F606W 



P(Ia|D)=0.05 
P(Ibc|D)=0.81 
P(II|D)=0.14 
z sn =0.886± 0.097 

Best-fit Model: 

Ic(SDSS-004012) 

MJD pk =55879 

z=0.886 

A v =0.00 

R v =3.1 

Am=0.3 

Xv=0-6 


EGSllPer 

<>F160Wil F125W# F350LP 



P(Ia|D)=0.99 
P(Ibc|D)=0.00 
P(II|D)=0.01 
z sn =0.915± 0.001 

Best-fit Model: 
Ia(SALT2) 
MJD pk =55702 
z=0.915 
Xl =-0.16 
C=0.1 

M-i 

xH-3 


COS12Mik 

Of160WD F125W^ F814W^ F606W 


P(Ia|D)=0.00 

P(Ibc|D)=0.13 

P(II|D)=0.87 



GND12Cal 

^F160WD F125W^ F814W# F350LP 


P(Ia|D)=1.00 

P(Ibc|D)=0.00 

P(II|D)=0.00 



GND13Wol 

♦ F160WD F125W- F814W# F350LP 



P(Ia|D)=0.31 
P(Ibc|D)=0.04 
P(II | D)=0. 65 
z SN =0.943± 0.001 

Best-fit Model: 

IIN(SDSS-012842) 

MJD pk =56496 

z=0.943 

A v =0.00 

R v =3.1 

Am=-0.2 

Xv=0-4 


GSDllMad 

<>F160Wn F125W^ F814W* F350LP 



P(Ia | D)=l .00 
P(Ibc | D)=0.00 
P(II|D)=0.00 
z SN =0.987± 0.001 

Best-fit Model: 
Ia(SALT2) 
MJD pk =55542 
z=0.987 
Xl =-1.42 
0 - 0.1 
3=4.1 

xH-3 


Figure 16. STARDUST light curve matches for the 16 SN from the CANDELS fields with redshifts 0.7 < z < 1.0, as in Figure 5. 
(A color version of this figure is available in the online journal.) 
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GNW12Gor 

OF160W F105W-F814W 



P(Ia|D)=0.09 

P(Ibc|D)=0.16 

P(II|D)=0.75 

z sn =1.016± 0.001 

Best-fit Model: 

IIN(SDSS-012842) 

MJD pk =56076 

z=1.016 

A v =0.48 

R v =3.1 

Am=-0.2 

xH-2 


GND13Reh 

Of160WD F125W* F350LP 



P(Ia|D)=0.00 
P(Ibc|D)=0.18 
P(II | D)=0. 82 

z sn =1.020± 0.001 

Best-fit Model: 
IIP(SDSS-017862) 
MJD pk =56421 
I z=1.019 
A v =0.00 
R v =3.1 
Am=-0.6 
XH-9 


GND13Jay 

#F160WO F125W^ F814W® F350LP 


P(Ia|D)=1.00 

P(Ibc|D)=0.00 

P(II|D)=0.00 



GND13Gar 

Of160WD F125W A F850LP# F350LP 



P(Ia | D)=l .00 
P(Ibc | D)=0.00 
P(II|D)=0.00 
z sn =1.067±0.0( 

Best-fit Model: 

Ia(SALT2) 

MJD pk =56424 

z=1.067 

xi=-1.42 

C=-0.1 

3=4.1 

Xv = 0.9 


COS12Her 

Of160WD F125W^ F814W- F606W 


P(Ia|D)=0.00 

P(Ibc|D)=0.07 

P(II|D)=0.93 



GSD12Agn 

O F 1 60W “ F850LP • F350LP 



P(Ia|D)=0.80 
P(Ibc|D)=0.02 
P(II | D)=0. 1 8 
z SN =1.095± 0.001 

Best-fit Model: 

Ia(SALT2) 

MJD pk =55966 

z=1.095 

Xl =-1.42 

C=0.3 

3=4.1 

xH-7 


EGS13Tan 

Of160WO F125W- F814W- F606W 



P(Ia|D)=0.15 
P(Ibc|D)=0.82 
P(II|D)=0.03 
z SN =1.137± 0.439 

Best-fit Model: 

Ic(SNLS-04D4jv) 

MJD pk =56415 

z=1.137 

A v =0.00 

R v =3.1 

Am=0.1 

xH-2 


GND12Tom 

Of160WO F125W^ F814W- 


P(Ia|D)=0.86 

P(Ibc|D)=0.14 

P(II|D)=0.00 



EGS13Wai 

^F160WO F125W^ F814W- F606W 


P(Ia|D)=1.00 

P(Ibc|D)=0.00 

P(II|D)=0.00 



GND13Cam 

OF160WU F125W* F350LP 
Of140W^F814W 




GSDllWas 

Ol'ieowQ F125W^ F814W® F350LP 


P(Ia|D)=1.00 

P(Ibc|D)=0.00 

P(II|D)=0.00 


56310 56340 56370 56400 56430 56460 



GND12Whe 

OF160W® F125W- F814W* F350LP 



P(Ia|D)=0.00 
P(Ibc|D)=0.00 
P(II | D)=l .00 
z SN =1.335± 0.005 

Best-fit Model: 

IIP(SDSS-020038) 

MJD pk =56134 

z=1.335 

A v =0.00 

R v =3.1 

Am=-1.4 

Xv = 0.6 


EGS13Cha 

Of160WD F125W^ F814W 



P(Ia|D)=0.36 
P(Ibc|D)=0.01 
P(II|D)=0.64 
z sn =1.406± 0.021 

Best-fit Model: 

IIN(SDSS-012842) 

MJD pk =56465 

z=1.406 

A v =0.50 

R v =3.1 

Am=-1.7 

X?=3-2 


EGSllOba 

Of160WD F125W- F814W- F606W 


MJD : 56430 56460 56490 56520 56550 56580 



P(Ia|D)=0.91 
P(Ibc|D)=0.09 
P(II|D)=0.00 
z sn =1.409± 0.001 

Best-fit Model: 
Ia(SALT2) 
MJD pk =55676 
z=l .409 
xi=-1.42 
0 - 0.1 
3=4.1 
xH-4 


EGSllNix 

♦ F160WD F125W^ F814W^ F606W 


P(Ia|D)=0.40 

P(Ibc|D)=0.13 

P(II|D)=0.47 



GND12Kin 

<>F160Wn F125W- F814W* F350LP 



P(Ia|D)=0.00 
P(Ibc | D)=0.44 
P(II|D)=0.56 
z SN = E460± 0.065 

Best-fit Model: 

IIN(SDSS-012842) 

MJD pk =56128 

z=1.603 

A v =0.50 

R v =3.1 

Am=-l.l 

xH-6 


GND13Fin 

Of160WD F125W- F814W# F350LP 


P(Ia|D)=0.03 

P(Ibc|D)=0.29 

P(II|D)=0.68 



Figure 17. STARDUST light curve matches for the 17 SN from the CANDELS fields with redshifts 1.0 < z < 1.5, as in Figure 5. 
(A color version of this figure is available in the online journal.) 
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Table 10 

Host Galaxies of 51 Supernovae with z < 1.5 


SN 

R.A. 

(J2000) 

Decl. 

(J2000) 

d 

(") 

d 

(kpc) a 

Morph. b 

SED C 

£host 

(±) 

z Reference d 

GSDIOTum 

03:32:17.871 

-27:50:59.48 

2.96 

57.9 

d 

A 

0.124 

0.001 

Le Fevre et al. (2004) 

COS12Cli 

10:00:16.060 

+02:12:37.38 

1.89 

37.8 

d 

SB 

0.187 

0.001 

Trump et al. (2009) 

GSD12Roc 

03:32:06.368 

-27:47:26.63 

0.00 

0.0 

i 

SB 

0.346 

0.003 

Keck+DEIMOS (B. Mobasher, S. Jha) 

EGSllTho 

14:19:31.685 

+52:51:56.05 

0.82 

16.8 

di 

A 

0.354 

0.001 

Keck+LRIS (A. Filippenko), Gemini+GMOS (S. Jha) 

COS12Ken 

10:00:35.978 

+02:15:25.81 

0.87 

17.6 

di 

SB 

0.373 

0.001 

Lilly et al. (2009) 

GND13She 

12:36:09.907 

+62:14:05.79 

0.30 

6.1 

sd 

A 

0.473 

0.001 

Wirth et al. (2004) 

GSDllTay 

03:32:54.502 

-27:47:03.67 

0.56 

11.5 

d 

A 

0.535 

0.001 

Le Fevre et al. (2004) 

GSDllFor 

03:32:14.320 

-27:47:13.15 

0.30 

6.3 

d 

A 

0.578 

0.001 

Mignoli et al. (2005) 

GND13Bid 

12:36:41.417 

+62:11:42.53 

0.72 

14.7 

i 

A 

0.585 

0.001 

Cohen et al. (2000) 

COS12Eis 

10:00:47.275 

+02:11:50.04 

0.78 

15.9 

sd 

SB 

0.605 

0.001 

Keck+DEIMOS (B. Mobasher) 

UDSllGar 

02:17:39.631 

-05:11:37.00 

0.14 

3.0 

sd 

SB 

0.651 

0.001 

Keck+DEIMOS (B. Mobasher) 

GSDllWor 

03:32:10.730 

-27:48:07.14 

0.12 

2.6 

d 

A 

0.654 

0.001 

Mignoli et al. (2005) 

GSDllRoo 

03:32:31.581 

-27:46:12.71 

0.20 

4.1 

s 

SB 

0.655 

0.001 

VLT+FORS2 (B. Leibundgut) 

GSWllJac 

03:32:32.310 

-27:54:20.46 

0.16 

3.3 

d 

A 

0.659 

0.001 

VLT+Xshooter (J. Hjorth) 

COS12Tru 

10:00:38.293 

+02:11:35.60 

0.92 

18.7 

d 

A 

0.665 

0.001 

Keck+DEIMOS (B. Mobasher) 

COS12Rea 

10:00:32.021 

+02:14:15.43 

1.72 

34.4 

d 

SB 

0.679 

0.001 

Lilly et al. (2009) 

GSDllBuc 

03:32:28.714 

-27:52:32.00 

1.14 

23.1 

s 

A 

0.679 

0.001 

Le Fevre et al. (2004) 

GSDllHar 

03:32:30.570 

-27:45:18.35 

1.75 

35.0 

d 

SB 

0.681 

0.001 

Le Fevre et al. (2004), Mignoli et al. (2005) 

COS 12 Aid 

10:00:15.212 

+02:17:30.84 

1.41 

28.5 

sd 

SB 

0.731 

0.001 

Keck+DEIMOS (B. Mobasher) 

GSDllLin 

03:32:29.799 

-27:49:19.26 

0.48 

9.9 

sd 

SB 

0.734 

0.001 

Keck+DEIMOS (A. Filippenko, S. Jha) 

GSDllAda 

03:32:19.785 

-27:54:09.15 

0.93 

18.9 

sd 

SB 

0.735 

0.001 

Balestra et al. (2010), Popesso et al. (2009) 

COS12Mon 

10:00:26.746 

+02:15:14.02 

0.31 

6.4 

s 

SB 

0.777 

0.001 

Lilly et al. (2009) 

GND13Ful 

12:36:19.227 

+62:15:12.76 

0.26 

5.3 

sd 

SB 

0.783 

0.001 

Wirth et al. (2004) 

GND12Bur 

12:36:32.464 

+62:15:32.75 

0.52 

10.7 

i 

SB 

0.69 

0.4 

phot-z (T. Dahlen) 

GND12Daw 










GND13Vin 

12:37:06.308 

+62:15:18.08 

0.43 

8.9 

i 

A 

0.840 

0.001 

Cohen et al. (2000) 

GNW12Ger 

12:37:16.965 

+62:20:38.34 

0.52 

10.7 

s 

A 

0.850 

0.005 

Keck+DEIMOS (A. Filippenko) 

GND13Roy 

12:36:46.232 

+62:15:27.46 

0.34 

6.9 

i 

SB 

0.851 

0.001 

Cooper et al. (2011) 

GSDllVan 

03:32:19.004 

-27:47:18.33 

0.61 

12.6 

d 

P 

0.74 

0.15 

phot-z (T. Dahlen) 

EGSllPer 

14:18:28.248 

+52:42:45.07 

1.03 

20.9 

d 

A 

0.915 

0.001 

Barro et al. (2011) 

COS12Mik 

10:00:31.733 

+02:26:13.92 

0.72 

14.6 

sd 

SB 

0.927 

0.001 

Keck+DEIMOS (B. Mobasher) 

GND12Cal 

12:37:10.476 

+62:15:47.40 

0.28 

5.8 

s 

SB 

0.941 

0.005 

HST+WFC3 (B. Weiner) 

GND13Wol 

12:36:58.933 

+62:18:10.04 

0.14 

3.0 

sd 

SB 

0.9431 

0.0002 

Barger et al. (2008) 

GSDllMad 

03:32:18.767 

-27:52:42.45 

0.44 

9.1 

d 

SB 

0.988 

0.001 

Keck+LRIS (A. Filippenko) 

GNW12Gor 

12:36:20.619 

+62:08:44.91 

0.62 

12.7 

u 

A 

1.0164 

0.0005 

Wirth et al. (2004), Treu et al. (2005) 

GND13Reh 

12:36:54.445 

+62:11:52.43 

0.06 

1.3 

s 

SB 

1.019 

0.001 

Barger et al. (2008) 

GND13Jay 

12:36:41.370 

+62:11:29.57 

0.46 

9.4 

d 

A 

1.03 

0.01 

HST+WFC3 (B. Weiner) 

GND13Gar 

12:36:40.813 

+62:11:14.34 

0.19 

3.9 

u 

SB 

1.07 

0.02 

phot-z (T. Dahlen) 

COS12Her 

10:00:47.438 

+02:15:17.35 

0.41 

8.4 

s 

SB 

1.1 

0.4 

Keck+DEIMOS (B. Mobasher) 

GSD12Agn 

03:32:25.861 

-27:50:19.81 

0.58 

11.8 

d 

SB 

1.095 

0.001 

Vanzella et al. (2008) 

EGS13Tan 

14:20:02.010 

+53:00:17.04 

1.05 

21.3 

sd 

A 

1.39 

0.8 

phot-z (T. Dahlen) 

GND12Tom 

12:36:42.636 

+62:18:20.74 

1.36 

27.4 

d 

A 

1.140 

0.001 

Barger et al. (2008) 

EGS13Wai 

14:20:28.495 

+53:04:58.72 

0.37 

7.6 

di 

SB 

1.1662 

0.0004 

Barro et al. (2011) 

GND13Cam 

12:37:07.357 

+62:10:26.94 

0.05 

0.9 

sd 

SB 

1.222 

0.002 

HST+WFC3 (B. Weiner) 

GSDllWas 

03:32:20.856 

-27:49:41.48 

0.00 

0.0 

d 

SB 

1.30 

0.05 

HST+WFC3 (A. Riess) 

GND12Whe 

12:36:24.206 

+62:17:39.92 

0.24 

5.0 

sd 

SB 

1.335 

0.001 

Keck+DEIMOS (C. Papovich) 

EGS13Cha 

14:20:24.200 

+53:03:37.50 

0.00 

0.0 

sd 

A 

1.40 

0.05 

Keck+MOSFIRE (M. Cooper) 

EGSllOba 

14:20:32.666 

+53:02:48.10 

0.08 

1.8 

di 

A 

1.409 

0.001 

Keck+LRIS (A. Filippenko) 

EGSllNix 

14:20:48.607 

+53:00:26.56 

0.10 

2.0 

u 

A 

1.73 

0.51 

phot-z (T. Dahlen) 

GND12Kin 

12:37:13.013 

+62:16:30.86 

0.06 

1.3 

u 

SB 

1.9 

0.4 

phot-z (T. Dahlen) 

GND13Fin 

12:37:18.322 

+62:17:09.34 

1.08 

22.0 

i 

SB 

1.4878 

0.0005 

Wirth et al. (2004), Reddy et al. (2006) 


Notes. 

a Physical separation between the SN and center of the host, computed from the measured angular separation in the preceding column, assuming a flat ACDM 
cosmology with Hq = 70, Q m = 0.3. 

b Visual classifications for host galaxy morphology: s = spheroid, d = disk, i = irregular. 
c Template-matching classification of host galaxy SED: P = Passive, A = Active, SB = Starburst type. 

d Unpublished spectroscopic observations are given as Observatory+Instrument ( name of PI). Host galaxy photometric redshifts are marked as phot-z (T. Dahlen et al., 
in preparation). 


two priors could be combined into a redshift-dependent, host- 
based prior — but that is beyond the scope of this paper. 

In Column 5 of Tables 3 and 9 we have reported the 
STARDUST probabilities derived using this host galaxy based 


prior. As one should expect, SNe for which we have abundant 
spectroscopic and photometric information are barely affected. 
Thus, classification probabilities that were close to 0 or 1 using 
the redshift-dependent prior typically do not shift. For objects 
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with intermediate probabilities, the la classification probability 
is almost uniformly increased, but the change is mostly within 
the range allowed by the classification uncertainties. The total 
change in the count of observed SNe la and the resulting 
change in the volumetric SN rate are presented in Table 8. 
Again, the shift is uniformly positive, but relatively small when 
compared to systematic uncertainties. The final column of this 
table reports the change in the SN la rate measurement as a 
fraction of the systematic uncertainty estimate, and we see that 
it is less than unity in every redshift bin. While not definitive, 
this result suggests that our baseline rates are not heavily biased 
by the redshift- dependent class prior, and any existing bias has 
been sufficiently accounted for in our systematic uncertainty 
estimates. 

APPENDIX B 

SUPPLEMENTARY DATA FOR SNe AT z < 1.5 

Discovery images for the 51 SNe with redshift z < 1.5 are 
shown in Figures 13 and 14. The names, positions, classification 
probabilities, and redshifts of those 51 low-z SNe are given in 
Table 9 (as in Table 3). Host galaxy information is provided in 
Table 10 (as in Table 4). Light curves and best-fitting template 
matches from STARDUST are shown in Figures 15-17. 
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